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: 1 Introduction 



The goal of these lecture notes is to introduce the developing research area of 
gravitational-wave phenomenology. In more concrete terms, they are meant to 
provide an overview of gravitational-wave sources and an introduction to the in- 
I terpretation of real gravitational wave detector data. They are, of course, limited 



in both regards. Either topic could be the subject of one or more books, and 



2 ■ certainly more than the few lectures possible in a summer school. Nevertheless, it 



is possible to talk about the problems of data analysis and give something of their 



Q I flavor, and do the same for gravitational wave sources that might be observed in 

I , the upcoming generation of sensitive detectors. These notes are an attempt to do 

5-1 . 

bJ)- just that. 

Despite an 83 year history, our best theory explaining the workings of gravity 
^ I — Einstein's theory of general relativity — is relatively untested compared to 

^ ' other physical theories. This owes principally to the fundamental weakness of the 

gravitational force: the precision measurements required to test the theory were 
not possible when Einstein first described it, or for many years thereafter. 

The direct detection of gravitational-waves is a central component of our first 
investigations into the dynamics of the weakest of the known fundamental forces: 
gravity. It is only in the last 35 years that general relativity has been put to 
significant test. Today, the first effects of static relativistic gravity beyond those 
described by Newton have been well-studied using precision measurements of the 
motion of the planets, their satellites and the principal asteroids. Dynamical 



* Supported by National Science Foundation awards PHY 98-00111 and PHY 95-03084 to The 
Pennsylvania State University, and PHY 93-08728 to Northwestern University, 
t Also Department of Physics, and Department of Astronomy and Astrophysics 



gravity has also been tested through detailed and comprehensive observations of 
the slow, secular decay of the Hulse- Taylor binary pulsar system.Bi What has not 
heretofore been possible is the direct observation of the effects of dynamical gravity 
on a laboratory instrument: i.e., the direct detection of gravitational radiation. 

The scientific importance of the direct detection of gravitational-waves does 
not stop with its detection, however. Strong gravitational- waves are difficult to 
generate: so difficult, in fact, that there is no possibility of a gravitational-wave 
"Hertz" -type experiment, where both the source and receiver are under laboratory 
control. The strongest gravitational-waves incident on Earth, as measured by our 
ability to detect them in the sensitive detectors now under construction, arise 
from astronomical sources. These are also the only sources that we can hope to 
observe in our detectors. The strongest of these anticipated sources — inspiraling 
or colliding neutron stars or black holes — are, in fact, of cosmological origin. 

Very little relevant detail is known about the gravitational-wave sources that 
we anticipate may be detectable in the instruments now under construction. Esti- 
mates of source strengths and event rates are difficult to make reliably. This is be- 
cause, at a deep and fundamental level, our understanding of the cosmos is limited 
to what we can learn from photons. The mechanism by which gravitational- waves 
are generated, on the other hand, favors sources that either do not radiate elec- 
tromagnetically {e.g., black holes), are obscured from view {e.g., the gravitational 
collapse of a stellar core), or are so distant and decoupled from the immediate 
origin of the corresponding electromagnetic radiation that we cannot reliably de- 
cipher the relevant source characteristics from the photons that reach us {e.g., 
7-ray bursts).! 

Gravitational-wave observations thus add a new dimension to our ability to 
observe the Universe: the observations that we make will tell us things we don't 
already know through other means. 

In order to describe sensibly the signature of gravitational-wave sources in real 
detectors we must first discuss in some detail how we characterize gravitational- 
waves, how we characterize gravitational-wave detectors, and how we give opera- 
tional meaning to the word "detect". These three subjects are addressed in sec- 
tions §1^ and respectively. In the context of gravitational-wave detection, 
gravitational-wave signals fall fairly neatly into three categories: burst signals, 
periodic signals and stochastic signals. Sources thought to be responsible for de- 
tectable signals in these categories are described in sections § [01 , § pl2| and § ^.3| , 



respectively. 



1.1 Conventions 

• The distance from detector to source will always be large compared to either 
a wavelength of the radiation field or the physical dimension of the detector; 
consequently, the incident radiation is effectively planar. 

• We choose a sign convention for the line element of Minkowskii spacetime 
and recall the Einstein summation convention, wherein repeated Greek in- 
dices in a product are implicitly summed over their full range: 

ds^ — ri^ijdx^dx^ (1) 

3 

= E V^udx'^dx'' (2) 

fj.,u=0 

= -c^ [dx°) ^ + [dx^) ^ + [dx^) ^ + {dx^) ^ (3) 
= -c^dt^ + dx^ + dy^ + dz^. (4) 

• We will always treat the gravitational fields as weak and use coordinates on 
spacetime that are either "Cartesian" -|-time or "spherical-polar" -|-time. 

• We will find it convenient to introduce the use of Latin indices to represent 
just the spatial components {i.e., x, y, and z) of a tenser. We generalize 
the Einstein summation convention to apply to repeated Latin indices in 
a product expression, with the implicit sum running over just the spatial 
coordinates, e.g., 

3 

x'yi = J2^^yr (5) 

• Except in the first several sections of these lecture notes, we will always work 
in units where the Newtonian gravitational constant G and the speed of light 
c are numerically equal to unity and the appearance of these constants in 
various formulae will be suppressed. Dimensional analysis will always suffice 
to determine how the formulae should appear with the constants in place. 

• Using c we can express time in units of length and frequency in units of in- 
verse length; similarly, by exploiting G and c we can express mass and energy 
in units of length. Power is then a dimensionless number. For CGS units. 



the conversion factors between mass, energy and length, and the physical 
constant with units of power, are 



G/c^ = 7.42 X 10"^^ cm/gm (6) 
G/c^ = 8.26 X 10-^° cm/erg (7) 
cVG = 3.63 X 10^9 erg/s. (8) 



2 Characterizing Gravitational Radiation 

For our purpose here — recognizing gravitational waves incident on a detector — 
two different characterizations of gravitational radiation are useful. The first is the 
radiation waveform and the second is the signal "power spectrum" . The waveform 
describes the radiation field's time dependence while the power spectrum describes 
its Fourier components. In § |2.1| and §2.3| we describe these different characteri- 
zations of gravitational radiation. Several important physical insights regarding 
gravitational radiation sources can be gained by considering the instantaneous 



power radiated by a source: we discuss these insights in §pl2 



2.1 Radiation waveform 

In this subsection we review briefly the expression of the radiation incident on a 
detector. Much of this section is by way of review; for more details, see either 
the lectures by Bob Wagoner in this collection, one of the many text books on 
relativity, or an excellent review article on the subject.! 

Gravitation manifests itself as spacetime curvature and gravitational waves as 
ripples in the curvature that appear to us, moving through time, to be propagating. 
Detectors are generally not directly sensitive to curvature, but to the mechanical 
displacement of their components; so, we focus our attention on the spacetime 
metric, from which physical distances between points in spacetime are determined. 
(The curvature is a function of the metric's second derivatives.) 

We assume that gravity is weak in and around our detector; correspondingly, 
we treat the spacetime metric as if it were the metric of Minkowskii spacetime, 
plus a small perturbation: 



where 7]^,^ is the Minkowski! metric and hf^„ the metric perturbation. The corre- 
sponding Une element, describing the proper distance between nearby spacetime 
events whose coordinate separation is the infinitesimal dx^, is 

ds^ = g^^dx^dx"" = ■q^^dx^dx'' + h^udx^dx"" . (10) 

Detecting gravitational waves amounts to building instruments that are sensitive 
to the effects of the small perturbation h^^; determining the signature of the 
gravitational waves in the detector thus requires determining h^i, and its influence 
on the detector. 

The metric g^j^i, tells us how the proper distance between points in spacetime 
is associated with our choice of coordinate system. Since the gravitational flelds 
near our detector are weak and the spacetime nearly Minkowskii, we can intro- 
duce coordinates that are, in the neighborhood of the detector, nearly the usual 
Minkowskii-Cartesian coordinates, with the deviations of the order of the pertur- 
bation. 

Now, small changes in the coordinate system do not change the proper dis- 
tance between events, only our labeling of them. If we make small changes in 
our coordinate system, of the order of the perturbation, then we will make cor- 
responding changes in the perturbation h^y. We can use this freedom to simplify 
the expression of h^^. Coordinate changes do not change the physics or any ob- 
servable constructed from /i^,^, of course. For this reason, and in analogy with 
electromagnetism, coordinate choices like these are referred to as gauge choices. 

With the separation of the metric into the Minkowskii metric 77^,^ plus a small 
perturbation, the fleld equations of general relativity become (at flrst order in the 
perturbation) a set of second order, linear differential equations for the ten com- 
ponents of the symmetric h^i,. Consequently, fixing the coordinates allows us to 
impose eight conditions on the ten components of the symmetric /i^i^, leaving just 
two dynamical degrees of freedom. These are identified as the two polarizations 
of the gravitational radiation field. 

An important gauge choice, always possible for radiative perturbations about 
Minkowskii space, is the Transverse-Traceless, or TT, gauge. Transverse- Traceless 
gauge is always associated with a particular observer of the radiation. Let the 
four- velocity of this observer have components U^. Without loss of generality let 
t mark the proper-time of this observer (so that Lf^ is just the coordinate vector 
in the t direction) and x, y and z be the usual Cartesian coordinates (to 0{h)) in 



the neighborhood of the observer. In TT-gauge the field equations are 

/ 02 \ 



subject to the constraints 



^ + V^1 V = -167rr^,, (11) 



h,M^ = 0, (12) 
^ {hjf.v'') = h,\k = 0, (13) 
h.kS'" = 0; (14) 



The metric perturbation satisfies a wave equation whose source is the stress- energy 
density T^^^ of the matter and (non-gravitational) fields. In more physical terms, 
the constraints are (in order): 

• hfj,^ is purely spatial; 

• the (spatial) metric perturbation hij is purely transverse: i.e., if the radia- 
tion wavevector is /c* (where the index i runs over just the spatial dimensions; 



see §[1.1D, then hijk^ vanishes for all j; and 



• the metric perturbation is trace-free. 

When there might be confusion we denote a metric perturbation in TT-gauge with 
a superscript TT on h^y] also, since the perturbation is purely spatial, we generally 
refer just to the spatial projection (in the coordinate system of the observer) of 
the perturbation, as in hj^ . 

Given a metric perturbation /i^^, expressed in any gauge, corresponding to 
a plane wave propagating in the direction n'^, we can recover the corresponding 
metric perturbation in TT-gauge by applying the linear operator Pi^ to the spatial 

UTT _ pU/ pm _\p plmr^l f^r\ 
Plm = Sim - niUm- (16) 

Here and henceforth we will always express the metric perturbation in TT-gauge. 

As mentioned above, gravitational wave detectors work by sensing the relative 
motion of their components induced by a passing gravitational wave. Let's see 
how the TT-gauge metric perturbation is related to such relative motion. 



Consider a single isolated test mass, initially at rest at coordinate position xa 
in a TT-gauge coordinate system. No forces act on this test mass; so, it moves 
through spacetime in such a way that its four-velocity always remains tangent 
to itself. (Forces, of course, cause the four- velocity to change direction.) The 
corresponding equations of motion for the spatial coordinates of the test mass are 



d f dxA \ _^ -pi dx^ dxj^ 
dr \ dr ) dr dr 



where r is the proper time of the test mass (initially r is equal to t since the test 
mass is at rest) and F^^ is the metric connection 

r^-y = Ig"^ (^M/3,7 + 9^^,0 - g^fs,,.) ■ (i8) 

(Recall that represents the derivative d/dx^.) Since the test mass is initially 
at coordinate rest the dx'^/dr vanish initially; so, the only connection component 
of interest is F^^. In TT-gauge, however, this component of the connection is 
identically zero (recall that /i^J is purely spatial); so, the equations of motion 
reduce to 

i.e., a free test particle at (TT-gauge) coordinate rest remains at coordinate rest. 
This apphes equally well for a second component of the detector, located at xb- 
it, too, remains at coordinate rest. 

This may seem paradoxical: if the coordinate separation of any two compo- 
nents of a detector remain unchanged by the passage of a gravitational wave, what 
is there to show the wave's existence? The paradox vanishes when we realize that 
coordinate separation is not physical separation. To determine the physical sep- 
aration of the detector's components we must invoke the metric again. Let the 
coordinate separation between the two components of the detector at time t be 
the infinitesimal dx^^j 

dx\s — x^Q — x\. (20) 

(Of course, dt — Q since we are talking about separation at the same coordinate 
time.) The physical distance between these two neighboring points in spacetime 
is 

ds^ = g^.^dx^'dx" (21) 
= rjuvdx^Bdx'^AB + dx^Asdx'AB (22) 
= rjjkdxABdx\B + hj^ dx^BdxAB- (23) 



The second term in equation ^ shows the effect of the gravitational wave on the 
separation between the two elements of the detector: as /i^i, oscillates, so does 
the distance. If the equilibrium separation between the components is L in the 
direction n-', to 0{h) the net change 5L in the separation is equal to 

5L = \LhJ^h^h'. (24) 

The physical distance between detector components does change, in an amount 
proportional to the undisturbed separation and the wave strength as projected on 
the separation. Gravitational wave detectors are designed to be sensitive to this 
displacement of their components. 

As mentioned above, the TT-gauge conditions amount to eight constraints on 
the ten otherwise independent components of the (symmetric) h^^,. There are thus 
two components of that are independent of the choice of coordinate system; 
correspondingly, in general relativity there are two dynamical degrees of freedom 
of the gravitational field. To see what what these amount to, without loss of 
generality consider a plane wave propagating in the z direction. Then we can 
write 

h iJJ dxf'dx'' = h+{x\t) (dx'^ - dy'^^ + 2h^{x^ ,t)dx dy, (25) 

where and /ix are the two independent dynamical degrees of freedom, or po- 
larizations, of the gravitational radiation field. 

Solutions to the wave equation for hij (eq. |lT]) can be analyzed in a slow motion 
expansion in exactly the same way as solutions to the Maxwell equations.0i'0 
The radiative hij in this expansion divide neatly into two classes of multipolar 
fields, which are (in analogy with electromagnetism) termed electric and mag- 
netic multipoles. The electric multipolar radiative fields are generated by time- 
varying multipole moments of the source matter density in the same way that the 
analogous electric moments of the Maxwell field are generated by the time- varying 
moments of the electric charge density. Similarly, the magnetic radiative moments 
are generated by the time-varying multipole moments of the matter momentum 
density, which is the analog of the electric current density. 

In electromagnetism, the first radiative moment of a charge distribution is a 
time- varying charge dipole moment. When electrical charge is replaced by gravi- 
tational charge — i.e., mass — we see that the corresponding dipole is just the po- 
sition of the system's center of mass, which (owing to momentum conservation) is 



unaccelerated. Consequently, in general relativity there is no gravitational dipole 
radiation. The first gravitationally radiative moment of a matter distribution 
arises from its "accelerating" quadrupole moment. Dotting the z's and crossing 
the f s, we find that, at leading order, the radiation field at a distant detector is 
related to the matter distribution of the source according to 



^^^-) = -^^.Q¥(t-r) (26) 

Qj^ = Pik{x)QkiPij{^)-\Pij{^)QimPim{x) (27) 

Qij{t) = J d^x p{t,x) (^XiXj -^6ij^ , (28) 
Pjk{x) = 6jk - XjXk/x'^. (29) 

The expression for hjj'" given above is the famous "quadrupole formula" of general 
relativity, which relates the acceleration of a source's quadrupole moment to the 
gravitational radiation emitted. It is, for weak gravitational fields, the exact 
analog of the more famous "dipole formula" of electromagnetism. 



2.2 Radiated power or energy 

Gravitational radiation carries energy away from the radiating system. Important 
insights into gravitational radiation can be gained by considering the energetics 
of radiation sources, which we do in this section. 

The instantaneous power carried by the radiation is, in the usual way, pro- 
portional to the square of the time derivative of the field integrated over a sphere 
surrounding the source: 

L oc ^Anr^h^ (30) 

The "exact" Q expression for the power carried away in electric quadrupole radia- 
tion is 

_ IG d'Qi.d'Q^A 
^-5?\^^/' ^^^^ 

where the <> indicates an average several periods of the radiation. Note that the 
power depends on Qij and not Qj/^- 

If we focus on the radiation emitted by a weak-field, dynamical source, we 
can use the multipolar expansions described above to replace the fields by the 



''In the context of our approximation of everywhere weak gravitational fields. 



multipole moments of the source. For an electric i-pole field radiated by a matter 
source with mass M, typical dimension R and internal velocity V, 



l(') oc - 
electric Q 



similarly, for a magnetic £-pole field 



GM fV\ 



cm \ cJ 



e+r 



(32) 



L^'^ oc - 

magnetic q 



GM fVV^^ 



cm 



(33) 



The total power radiated is the sum over the power radiated in each of the mul- 
tipoles. 

Aside from numerical factors and symmetries, power radiated in the electric 
^-pole channel is suppressed relative to that in the electric quadrupole channel by 
a factor of {V/CY'^^^'^\ Similarly, the radiation in the magnetic £-pole channel is 
suppressed from the electric quadrupole radiation by a factor of {V / GY^^^^\ Con- 
sequently, sources whose internal velocities are significantly less than the speed 
of light radiate principally in the electric quadrupole channel (again, unless sup- 
pressed by symmetries). 

There is still another way of looking at the power radiated by a gravitational 
radiation source. For the gravitational wave detectors we can hope to build all the 
radiation of interest is of astrophysical origin. Excepting only a stochastic grav- 
itational wave background, the radiation sources are all distinct systems whose 
structure or dynamics are governed by gravity. For these systems, judicious ap- 
plication of the Virial TheoremEl allows us to relate the internal velocities V to 
the depth of the gravitational potential GM/R, 



- . (34) 



Thus, for astrophysical sources 



e+3 



^magnetic G U 7 ~g{^RJ ' ^ ' 

Strong gravitational wave sources thus have strong internal gravitational fields. 

Finally, dimensional analysis of equation ^ for the power radiated in the 
electric quadrupole leads to an important physical insight. Dimensionally, the 



system's quadrupole moment is proportional MR^ . In a closed, radiating system 
there is a typical time scale T for motion within the source; consequently, the 
total radiated power can be written 

MR^/T^ MR^ , , 

The quantity MR'^/T'^ ^ MV"^ can be interpreted as the kinetic energy of source 
matter engaged in motion associated with a time-varying quadrupolar moment. 
Similarly, we identify MR'^/T^ as the instantaneous power available to be radiated. 



Not all this power is radiated, however. Equation ^ shows that the fraction of 
the available power actually radiated is equal to the available power divided by a 
"fundamental power" defined by the physical constants G and c: 

— = 3.6 X lO^^erg/s. (38) 

Gr 

The magnitude of this fundamental power gives us a feeling for the weakness 
of the gravitational interaction. For a source to radiate even one part in 10^ of 
the power available to be tapped by the radiation field, it must have internal 
motions where the kinetic energy involved in quadrupolar motion is greater than 
3.6 X 10^° erg/s. For scale, this is 10^"^ times greater than the power liberated in 
all of the nuclear reactions occurring in the Sun! 



2.3 Signal Power Spectrum 

Observations of gravitational wave signals are always of finite duration: either the 
signal is a burst of duration less than the observation period or the signal duration 
is determined by the period between when the detector is turned on and when 
it is turned off. A useful characterization of this observed signal is its spectrum: 
the contribution to the overall mean square signal amplitude owing to its Fourier 
components at a given frequency. 

For definiteness, focus attention upon some particular polarization h{t) of a 
gravitational-wave signal that is observed over a period T beginning at t = 0. The 
Fourier transform of this signal is 

h{f) = r dte-^^'f%{t). (39) 

JO 

tWe use the engineering convention for the Fourier transform. 



The signal spectrum is evaluated for positive frequencies and is twice the square 
modulus of its Fourier transform averaged over the observation, or 

W):=|lM/)r (40) 
for non-negative /. Since h{t) is real, we can use The Parseval Theorem to obtain 

I dfPhif) = fl dt\h{t)\' = {h'), (41) 

where (■) denotes a time average. The signal spectral density is thus a measure of 
the contribution to the mean-square signal amplitude owing to its Fourier com- 
ponents in a unit bandwidth. (For non-burst — i.e., stochastic or periodic — 
signals, we often take the limit as T ^ oo.) 

As we have described it, the signal spectrum is derived from the signal wave- 
form h{t) by "throwing away" the phase information. There is clearly much less 
information in P{f) than in the corresponding h{t): why, then, is P{f) an inter- 
esting characterization of a signal? 

One reason is that real detectors are only sensitive to radiation in a lim- 
ited bandwidth — i.e., at certain frequencies. The integral of the signal power 
spectrum over the detector bandwidth is the contribution to the mean-square 
amplitude of h from power in the detector bandwidth. 

A second reason is that it is not always possible to determine the waveform 
of a gravitational wave signal. For example, the waveform of a stochastic signal, 
arising from a primordial background or from the confusion limit of a large number 
of weak sources, is intrinsically unknowable. Nevertheless, the signal spectrum is 
straightforward to calculate. In this case, the spectrum embodies everything we 
can know about the gravitational wave signal. 

Another example illustrates a different circumstance. Calculations of gravita- 
tional radiation waveforms /i+(t) and hy^{t) from the kind of stellar core collapse 
that triggers type II supernovae are, even in their grossest details, extremely sensi- 
tive to the details of the stellar model and the physics included in the simulations. 
In the face of this variety of structure, however, the spectra all show a remark- 
able similarity.0 It may be that this variety reflects our ignorance of the relevant 
physics and that with better understanding the waveforms would show much less 
variation and much greater predictability; it may also be that the details of the 
collapse waveform are in fact very sensitive to the initial conditions. Whether in 



practice or in principle, the waveform is today unknown; nevertheless, the spec- 
trum does appear to characterize the signal quite well. 

We close with a final reason that the spectrum is a useful characterization of a 
gravitational wave signal. The sensitivity of a gravitational wave detector is lim- 
ited by the detector noise, which is an intrinsically stochastic process. In the best 
detectors, the noise is fully characterized by its spectrum (cf. |3.5| ). We expect 
intuitively that a signal is detectable only when its spectrum has greater magni- 
tude than the the detector noise spectrum over a sufficient range of frequencies. 
This qualitative notion finds quantitative expression in the signal-to-noise ratio. 



which we discuss in 53.5 below. 



2.4 Conclusion 

For the purposes of detection, gravitational waves are usefully characterized by 
their waveform or spectrum. There are important sources for which the explicit 
waveform is not known, either because it is intrinsically unknowable, our grasp of 
the underlying physics is not complete or the calculations involved in determining 
it our beyond our capabilities. In these cases, it may still be possible to estimate 
the signal spectrum, which then serves to characterize it. 



3 Characterizing The Detector 
3.1 Introduction 

Gravitational- wave detectors transform incident gravitational waves into, e.g., 
electrical signals that we can more easily manipulate. In § |3.2| , we describe briefiy 
and schematically two of the detector technologies currently being pursued to de- 
tect gravitational waves. For all detectors we might realistically imagine building 
the detector response is linear in the incident radiation: i.e., the time history of 
the detector output is linearly related to the time history of the incident radiation. 
There are two aspects of this response that we must consider: differential sensitiv- 
ity to the radiation incident from different directions, and differential sensitivity 
to incident radiation of different frequencies. The first of these is described by the 
detector's antenna pattern, which we discuss in §^.3|, and the second of these is 



described by the detector's response function, which we discuss briefiy in §^75 



The output of a gravitational wave detector might contain a particular grav- 
itational wave signal; however, it always contains noise. Detection, discussed in 
§^ below, involves distinguishing gravitational wave signals from detector noise. 
To make this distinction we must have some characterization of the signal (e.g., 
by waveform or by spectral density) and detector noise. How we characterize 
detector noise is the subject of § |3.5| . 



3.2 Gravitational Wave Detectors 
3.2.1 Acoustic Detectors 

The earliest and most mature detection technology is, conceptually, nothing more 
than a high quality tuning fork. Gravitational waves excite the tuning fork; grav- 
itational waves at or near the tuning fork resonant frequency excite it into large 
amplitude oscillations. The tuning fork is instrumented so that its acoustic vibra- 
tions become electrical signals, which, when amplified, are the gravitational wave 
signal. 

Physically, the tuning fork is realized as one or more normal modes of a large 
metal-alloy block: the fundamental longitudinal mode of a right cylinder for the 
currently operating detectors, the five quadrupole modes of a sphere or a truncated 
icosahedron0 for the proposed next generation detectors. The choices made in the 
construction of the ALLEGRoSlli detector, built and operated at the Louisiana 
State University, are typical for the current generation of these right cylindrical 
"bars": diameter of 60 cm, length of 3 m, and cast of A15056 alloy for a total 
mass of 2296 Kg. 

The mechanical oscillations of the tuning fork are converted into electrical sig- 
nals, which are then amplified, digitized, and otherwise manipulated to determine 
whether gravitational waves are present or absent. In all of the high-sensitivity 
bar detectors operating today, the transducer is not directly connected to the bar, 
but instruments a second mechanical oscillator, of lower mass and smaller physical 
dimension, that is itself coupled to the bar. Gravitational waves drive the bar, 
which in turn drives the second oscillator. In the process, the amplitude of the me- 
chanical vibrations are amplified, and it is this mechanically amplified motion that 
is converted into electrical signals and further amplified, etc. The coupling of the 
two mechanical oscillators splits the fundamental longitudinal mode of the larger 
bar into two closely-spaced modes. For the ALLEGRO detector, the antenna's 



normal modes are at 896.8 Hz and 920.3 Hz. 

At this writing there are five operating cryogenic acoustic gravitational wave 
detectors: 

• ALLEGRO, at the Louisiana State University in the United States,^ 

• AURIGA, at the University of Padua in Italy, 

• EXPLORER, operated by the Rome group and located at CERN,0 

• NAUTILUS, operated by the Rome group and located at the Frascati INFN 
Laboratory ,0 and 

• NIOBE, at the University of Western Australia.0 

In addition to these classical "bar" detectors, several spherical or truncated icosa- 
hedral detectors have been proposed or are undergoing technical development: 
SFERA, TIGA,0ffl GRAIL,il and OMNI.i 

3.2.2 Interferometric Gravitational Wave Detectors 

An alternative technology for the detection of gravitational waves involves the use 
of a right-angle Michelson interferometer with freely suspended mirrors. Gravita- 
tional waves incident normal to the plane of an interferometer will lead to differ- 
ential changes in the distance between the corner and end mirrors. For frequencies 
much less than the light storage time in an interferometer arm, the corresponding 
motion of the fringes is proportional to the incident radiation waveform. 

There are currently two Km-scale interferometer projects under construction: 
the French/Italian VIRGO Projecti and the US LIGO Project.^ VIRGO will 
consist of a single interferometer with 3 Km long arms situated just outside of Pisa, 
Italy. LIGO will consist of two separate facilities, one at Hanford, Washington 
and one in Livingston, Louisiana. Each LIGO facility will house an interferometer 
with arms of length 4 Km; in addition, the Hanford facility will also hold an 
interferometer of 2 Km arm length in the same vacuum system. 

In addition to these larger interferometer projects, there are three more in- 
terferometric detectors of somewhat smaller scale under construction: the Aus- 
tralian ACIGA project, the German/U.K. GEO 600 projectil and the Japanese 
TAMA 300 project. The ACIGA Project's ultimate goal is a multi-kilometer de- 
tector, to be located several hours outside of Perth; presently, they are beginning 



the construction of an approximately 80 m prototype at the same site. GEO 600, 
located in Hanover, Germany, is a folded Michelson interferometer with an op- 
tical arm length of 1.2 Km. The Japanese TAMA 300 is a 300 m Fabrey-Perot 
interferometer located just outside of Tokyo; it is hoped that the success of this 
project will lead to the construction of the proposed Laser Gravitational Radiation 
Telescope (LGRT), which would be located near the Super-K neutrino detector. 

There are several ways to make an interferometer more sensitive at frequencies 
less than the reciprocal of the detector's light storage time. One is to increase its 
arm length (recall equation |2j!). The Laser Interferometer Space Antenna — LISA 
— is an ambitious project to place in solar orbit a constellation of satellites that 
will act as an interferometric gravitational wave detector.@'0 The arm length of 
this interferometer would be 5 x 10^ Km. The LISA project has been approved by 
the European Space Agency as part of its Horizon 2000+ Program; additionally, 
the US National Aeronautics and Space Administration is actively considering 
joining ESA as a partner to accelerate the development and launch of this exciting 
project. 

3.3 Antenna Pattern 

Gravitational wave detectors respond linearly to the applied field. The interfero- 
metric and bar gravitational wave detectors now under construction or in opera- 
tion have only a single "gravitational wave" output channel.^ When a plane wave 
is incident on such a detector, the time history of the output channel is linearly 
related to a superposition h{t) of the + and x polarizations of the incident plane 
wave: 

h = F+h+ + F^h^. (42) 

The factors and Fx describe the detector's "antenna pattern" , or differential 
sensitivity to radiation of different polarizations incident from different directions. 
(In fact, the antenna pattern may also be a function of radiation wavelength; 
however, when the wavelength is much larger than the detector this dependence is 
insignificant.) They depend on relative orientation of the plane- waves propagation 
direction and the definition of the + and x polarizations 

•f Some proposed acoustic detectors are instrumented on several independent modes. In this case, 
each mode may be considered a separate detector and represented as a single gravitational wave 
channel. 



If we fix the propagation direction and rotate the polarization of the incident 
radiation, then the detector response h(t) will change. Define the polarization 
averaged root-mean-square (RMS) antenna pattern F, 

F\k) = Fl{k) + Fl{k), (43) 

where k is the wave-vector of the incident plane wave and the overline denotes 
an average over a rotation of the incident radiation's polarization plane. The re- 
sult depends only on the wave-vector (or, alternatively, the wave's propagation 
direction and wavelength) and is proportional to the detector's root-mean-square 
response to plane-wave radiation incident from a fixed direction at fixed wave- 
length. For either the acoustic or interferometric detectors now operating or under 
constructions, F{k) is independent of the magnitude of k as long as the radiation 
wavelength is much larger than the detector. 

A convenient pictorial representation of the detector's response results if we 
plot the surface defined by nF{k) for fixed where h is the unit vector in the 
direction of the source relative to the detector {i.e., n = —k/k). In such a figure, 
the response of the detector to a plane wave with wave-vector k (appropriately 
averaged over polarization) is proportional to the distance of the surface from 
the origin in the direction of the source (n). In the remainder of this subsection 
we describe the antenna pattern of interferometric and acoustic bar detectors to 
incident gravitational plane waves. 

3.3.1 Bar detectors. 

In a classic bar detector, incident gravitational waves drive the fundamental lon- 
gitudinal mode of a right cylindrical bar. The driving force — and thus the 
radiation — is determined by observing the motion of this mode. For definite- 
ness, let the longitudinal axis of the bar be along the i-direction, and consider a 
plane gravitational wave propagating in the x-direction: 

hj^{t)dx'dx^ = h+{t) (dy^ - dz^) + 2h^{t)dydz. (44) 

The -|- polarization mode changes the z-distance between the atoms in the bar. 
This change is resisted by inter-atomic forces in the bar; thus, the bar's longitu- 
dinal normal mode is driven by this polarization component of the incident wave. 
The X polarization, on the other hand, does not excite the bar's mode in this 
way; so, the detector is insensitive to this component of the incident radiation. 



Fig. 1. The polarization-averaged RMS sensitivity of a bar detector to gravita- 
tional waves incident from any direction. The detector is at the origin of the figure 
and has its symmetry axis along the figure's z axis. The magnitude of the dis- 
tance from the origin to the surface in a direction h is proportional to the relative 
response of the detector to radiation incident on the detector from that direction. 

If, on the other hand, the waves are incident along the longitudinal axis, then 
neither the x nor the -|- polarization components cause any change in the longitu- 
dinal distance between the atoms in the bar; correspondingly, the bar is insensitive 
to waves incident along on the bar along its axis. 

Finding the response to radiation of different polarizations incident from direc- 
tions intermediate between these extremes is a relatively straightforward exercise 
in geometry. First define the bar's coordinate system. Let the bar's symmetry axis 
define the z direction and choose x and y such that (x, z) defines a right-handed 
coordinate system. 

Next consider a plane gravitational wave propagating in the k direction and 
define the polarizations of the gravitational wave. Introduce a plane orthogonal 
to /c. In this plane, spanned by the Cartesian coordinates x' and y', the TT- gauge 



metric perturbation can be written 

hijdx'^dx-' = (^dx''^ — dy'"^^ + 2hxdx'dy' . (45) 

The choice of x' and y' is arbitrary: different choices correspond to either or both 
a rotation of one polarization state into another or a reflection that flips the sign 
of hx- For definiteness, we choose x' to be orthogonal to z and y' such that {x', 
y', k) is a right-handed coordinate system. (In the degenerate case, where k is 
parallel to i, we choose x' parallel to x and y' parallel to y.) With these choices, 
it is straightforward to show that the response of the bar's longitudinal mode is 
proportional to /i, where 

h = F+h+ + Fxhx (46) 
F+ = sin^^ (47) 
Fx = sin^^ (48) 
cos^^ = (k-hy (49) 

Figure ^ shows the polarization-averaged RMS sensitivity of a right cylindrical 
acoustic detector to plane waves incident from a given direction. To interpret the 
figure, imagine a detector at the figure's origin with its symmetry axis coincident 
with the figure's z axis. A plane wave, arriving from direction n, leads to a 
detector response proportional to the distance in the direction n from the figure 's 
origin to the surface. 

3.3.2 Interferometric detectors. 

Interferometric gravitational wave detectors respond when incident gravitational 
waves cause a differential change in the length of the interferometer arms. Focus 
attention on interferometers whose arms meet in a right angle. To get a sense of 
the differential sensitivity of such a detector to radiation of different polarizations 
incident from different directions, define a right-handed interferometer coordinate 
system whose origin is the intersection of the arms and whose x and y coordinate 
directions are in the direction of the arms. Let a plane wave, described by the 
perturbation 

hijdx^dx-^ = /i+ (dx^ — dy"^^ + 2hxdxdy, (50) 

be incident on the detector from direction z. There will be no detector output 
proportional to hx, since that component of the radiation does not lead to a 



differential change in the arm lengths; on the other hand, the polarization com- 
ponent proportional to hj^ does lead to a differential change in the arm lengths 
and, correspondingly, to detector output. 

Similarly, consider radiation incident on the detector along the interferometer's 
X arm: 

hijdx^dx^ = h+ {dy^ — dz^^ + 2hxdydz. (51) 

Again, the x polarization mode does not lead to a differential change in the 
interferometer arm lengths (at first order in h); so, the detector is not sensitive 
to radiation with this polarization. On the other hand, radiation in the + 
polarization mode, as we have defined it, leads to changes in the length of the y 
arm while leaving the x arm length unchanged; consequently, there is a differential 
change in the interferometer arm length and the detector is sensitive to radiation 
of this polarization incident from this direction. 

To determine in general the coefficients F+ and Fx that describe the response 
of an interferometric detector to incident plane waves, first describe the polariza- 
tion modes of radiation incident on the detector relative to the detector coordinate 
system. In the usual {9, 0) spherical coordinates associated with the interferom- 
eter coordinate system, the incident direction of a plane-wave propagating with 

— * 

wave-vector k is 

cos^ = -k-z/\k\, (52) 

tan0 = (53) 
k ■ X 

In the plane orthogonal to the radiation propagation direction k, let the x' direc- 
tion be parallel to the xy-plane and the y' direction be orthogonal to x' so that 
{x',y',—k) forms a right-handed coordinate system. [In the degenerate case — 
radiation propagating parallel to the z direction — we take x' parallel to x and y' 
such that {x'.y', —k) is right-handed.] In terms of this coordinate system, define 
the -|- and x polarizations of an incident gravitational wave by 

hijdx'dx^ = h+ (dx'^ - dy'^) + 2Kdx'dy'] (54) 

then, the antenna pattern factors F+ and Fx are given by 

F+ = ^ (^1 -|-cos^^^ cos20cos2'0 — cos^sin20sin2'0, (55) 

Fx = - (l -|- cos^^] cos20sin2'0 -|- cos6'sin20cos2'0. (56) 



Fig. 2. The polarization-averaged RMS sensitivity of an interferometric gravita- 
tional wave detector to radiation incident from any direction. The detector is at 
the origin of the figure and has its arms aligned with the figure's x and y axes. 
The magnitude of the distance from the origin to the surface in a direction n is 
proportional to the relative response of the detector to radiation incident on the 
detector from that direction, averaged over all polarizations. 

Figure |^ shows the polarization-averaged RMS sensitivity of a right-angle inter- 
ferometric detector to plane waves incident from a given direction. The detector 
is at the origin of the figure, with its arms along the figure's x and y axes. The 
detector's sensitivity to radiation incident on the detector from direction n is pro- 
portional to the distance of the surface from the figure's origin in the direction 
n. 

3.4 Response Function 

The output of a gravitational wave detector is a voltage, v{t), that is linearly 
related to the incident radiation. Consider a gravitational plane wave, with po- 



larizations /i+ and hx-, incident on a detector with antenna pattern described by 
and Fx. The detector response is given by 

v{t) = f dTh{T)K{t-T), where (57) 
h{t) = F+h+{t) + FxKit) (58) 

and K is the kernel of the hnear transformation. 

It is instructive to express this convolution in the frequency domain: 



v{t)= dfh{f)K*{f)exp{2mft), (59) 

J —oo 

where 

roo 

gif)= dtgit)exp{-2mft) (60) 



and we have assumed that K(t) vanishes for negative r.^ From equation ^ we 
see that the response of the system — the output voltage — depends on the 
frequency of the incident radiation: depending on the character of the detector, 
the response may be relatively large for some frequencies and relatively small for 
others. 

As an example, consider two equal masses M connected by a spring (spring 
constant uj^M, quality factor Q). Denote the equilibrium separation by L. A 
passing gravitational wave of appropriate polarization disturbs the equilibrium 
separation of the system. The net result is that the passage of a gravitational 
wave acts as a driving force on the system's normal mode: 

ot + -^x + UqX = -Lfi, (61) 

where x is the difference between the actual and equilibrium separation. 

Suppose that we instrument this system with strain gauges to produce an out- 
put voltage v{t) proportional to x{t), the deviation from equilibrium separation. 
How is v{t) [i.e., x{t)) related to h{t)'! In the frequency domain we see that 

where 

/o = ^o/27r. (63) 
§ Corresponding to a causal impulse response! 



The output voltage for excitations near the resonance can vary dramatically as a 
function of frequency. 

The response function we have just described is equivalent conceptually to that 
of a modern acoustic detector: the radiation manifests itself as a driving force on 
the system's normal modes and the response is a strong function of the frequency 
in the neighborhood of the resonances. 

The response function of an interferometric detector is quite different. For an 
interferometer, at frequencies much below the round-trip travel time (but greater 
than the pendulum frequency of the suspended mirrors and beam-splitter) the 
detector response is independent of frequency; only when the frequency becomes 
comparable to or larger than the round-trip light travel time in an interferometer 
arm does the response vary with frequency.^^ 

The amplitude of the response determines those frequencies where an incident 
gravitational wave of unit amplitude gives relatively large amplitude output and 
where it gives relatively small amplitude output. It is not, however, the case 
that relativity large amplitude output corresponds to relatively large sensitivity, 
if by sensitivity we mean greater ability to detect. To address the question of 
sensitivity, we must turn to yet a different aspect of a detector's function: its 
noise. 

3.5 Noise 

The output channel of a gravitational wave detector is always alive with ran- 
dom fluctuations — noise — even in the absence of a gravitational wave signal. 
In a perfect world noise would arise exclusively from fundamental physical pro- 
cesses: e.g., fluctuations owing to the finite temperature of the detector, counting 

lit is commonly said that an interferometer responds to a passing gravitational wave propor- 
tionately with the differential change in the IFO arm length. This is not quite right. The 
response of an interferometric detector to a passing gravitational wave is proportional to the dif- 
ferential change in the round-trip light travel time in the arms. The round-trip light travel time 
involves the integrated change in the arm length over the past, as opposed to the instantaneous 
separation at the time of reflections. For frequencies small compared to the inverse round-trip 
light travel time the difference is negligible. 

It is also the case for interferometers that the frequency dependence of the response function 
varies with the incidence direction of the radiation though — again — this is only significant at 
frequencies comparable to or larger than the inverse round-trip light travel time in an arm. 



statistics of individual photons on a photo-dctcctor, etc. In the less than perfect 
world in which we live there will be other contributions to the detector noise, 
beyond these fundamental processes, that arise from the imperfect construction 
of the detector {e.g., bad electrical contacts), imperfections in the materials used 
to construct the detectors {e.g., mechanical creep and strain release), and from 
the detector's interaction with the (non-gravitational wave) environment {e.g., 
seismic vibrations, electromagnetic interactions, etc.). 

Detection of gravitational waves requires that we be able to distinguish, in the 
detector output, between signal and noise. This requires that we have character- 
ized the noise (and not only the signal). Since noise is intrinsically random in 
character, that characterization is in terms of its statistical properties. Some of 
these statistical properties we can predict, model or anticipate a priori, based on 
the detector design; nevertheless, it is important to realize that an experimental 
apparatus is a real thing made in the real world and will never behave ideally. 
While a large part of the experimental craft involves building instruments that 
operate as close as possible to their theoretical limits or prior expectations, the 
final characterization of a detector will always be determined or verified empir- 
ically. In this section we describe something of how noise in gravitational wave 
detectors is characterized. 

3.5.1 Correlations 

Just as a probability distribution is fully characterized by its moments, so the 
random output of a gravitational wave detector can be fully characterized by its 
correlations. The A'"-point correlation function describes the mean value of the 
product of the detector output sampled at N different times. Mean, in this case, 
refers to an ensemble average, where the ensemble is an infinite number of identi- 
cally constructed detectors. Denoting by n{t) the noisy output of a gravitational 
wave detector in the absence of any signal, the A'"-point correlation function of the 
noise distribution is given by 

C7v(ro, . . . , tn-i) = n{To) . . . n(Tiv_i), (64) 

where the over-bar signifies an ensemble average, which is also referred to as an 
average across the process. 

As a practical matter ensemble averages are impossible to realize experimen- 
tally: one rarely has the opportunity of working with even two similar detectors. 



let alone an infinite number of identical ones. Thus, while a handy theoretical 
construct, the general set of correlation functions is not of great practical use in 
characterizing the behavior of a real detector. 



3.5.2 Stationarity 

If, however, the behavior of the detector noise does not depend significantly on 
time — i.e., the noise is stationary — then the utility of the correlation function 
as a practical tool for characterizing detector noise increases dramatically. When 
the noise character is, figuratively, the same today as it was yesterday and as 
it will be tomorrow, then the detector yesterday (or an hour, or a minute, or a 
second ago) can be regarded as an identical copy of the detector we are looking at 
now, and both are identical copies of the detector tomorrow. Consequently, in the 
spirit of the ergodic theorem, we can replace the average across the process — the 
ensemble average — with an average along the process — a time average. The 
iV-point correlation function is then a function of the difference in time between 
the samples: 

1 /•*o 

CAr(ri, . . . , tn-i) = lim — / n(t)n(t - n) ■ ■ ■ n(t - tn-i) dt. (65) 

Of course, perfect stationarity is an impossible requirement. As a practical 
matter, what we require is that the noise process be stationary over a suitably 
long period. Let's try to make that concept more quantitative. To simplify the 
discussion, assume (without loss of generality) that the noise process has zero 
mean. Consider first the two-time correlation function of a stationary process: 

C2(r) = lim 4 r n(t)n(t - r) dt. (66) 

T^oo T Jto-T 

For sufficiently large r we expect intuitively that 6*2(7") should vanish: the output 
now should be effectively uncorrelated with the output in either the distant past 
or the distant future. This will also be the case for the higher-order moments 
as well: for sufficiently large (any k), the correlation function Cn should van- 
ish. Thus, we don't need to require perfect stationarity; rather, we require only 
that the statistical character be approximately stationary, varying significantly 
only over times long compared to the longest correlation time. In that case, we 



can approximate the correlations Cat by averaging, as in equation |6^, over finite 
periods. 



3.5.3 Gaussian Noise 



Noise from fundamental processes tends to be either Gaussian [i.e., originating 
from contact with a heat bath or some dissipative process) or Poissonian {e.g., 
originating in the counting statistics of identical and independently distributed 
— i.i.d. — events that occur at a fixed, average rate). For the gravitational wave 
detectors under construction, the intrinsically Poissonian processes {e.g., photon 
counting statistics) have rates so high that they can be treated as Gaussian and 
we do so here and below. 

One way to think about the detector noise is as a superposition of a Gaus- 
sian and approximately stationary component, a (hopefully lower amplitude) non- 
Gaussian, but still stationary component, superposed finally with a non- stationary 
component. General statements cannot be made about the non-Gaussian or non- 
stationary components: they differ from instrument to instrument and environ- 
ment to environment and can only be characterized empirically. The characteriza- 
tion of the Gaussian-stationary component, however, is remarkably simple and has 
a useful physical interpretation, which we review in this and the next subsection. 
(For more information and detail, see Finn.0) 

Up to now we have considered the output of a detector as an analog process: 
i.e., one that is continuous in time. In fact, the output we observe will have been 
sampled discretely at some sampling rate fs, chosen to be something more than 
twice as great as the maximum frequency of interest for the detector output. So, 
instead of writing the noise at the detector output as n{t) we write 



for constant At. 

When the noise in the detector is Gaussian and stationary, any single sample 
n[j] of the detector output is drawn from a normal distribution with a mean and 
variance that are independent of when the sample was taken. Without loss of 
generality we can assume that the mean n vanishes, in which case 




(67) 



where 



tf^ = to + kAt 



(68) 



We understand the variance cr^ of the distribution to be the ensemble average of 
the square of the detector noise: 



a' = n\jf. (70) 

Equation ^ holds true for each sample n[j]; consequently, the joint probability 
that the length sequence of samples n[j], j running from 1 to N^, is a sample 
of detector noise is given by the multivariate Gaussian distribution 



exp 



-lJ:rk=on[j]\\C-%,n[k] 



P(n[l], n[2], . . . , n[Nr]) = (71) 



(27r)^^det||C| 

In place of the variance o"^ that appears in the exponent of equation ^ is the co- 
variance matrix C. (The matrix C^^ that appears in equation |7l| is, by construc- 



tion, positive definite; consequently, it is non-singular and invertible.) Similarly, 
in place of the factor that appears in the denominator of equation ^ is the 
determinant det 1 1 C 1 1 . 

The mean over the product is the value of the correlation function 

C2(tj — tk); it is also just the value of the jk element of the covariance matrix 

C,[j-k]=n[j]n[k] = \\C\\^,. (72) 

Since the detector noise is also assumed to be stationary, ||C||j.j;. can depend only 
on the difference j — k; correspondingly, C is constant on its diagonals: i.e., it is 
a Toeplitz matrix. Consequently, it is fully characterized by the sequence c[k] of 
length 2Nj- — 1 whose elements are the first row and column of C: 

c[j-k] = \ \C\\^, = C,{t,-tk). (73) 

The sequence c and the process mean (which we have assumed to vanish) fully 
characterize the random process. The sequence c, however, is just the two-time 
correlation function of the detector output! Thus, the two-time correlation func- 
tion C2 fully characterizes a Gaussian stationary process: all the higher order 
correlation functions Cat either vanish (for odd A^) or are expressible as sums of 
products of 6*2. Once we have determined C2, then, we have completely deter- 
mined the character of the Gaussian noise process. 



3.5.4 Likelihood function 



In the last section we evaluated 

^ probability of observing 

P(t;|0) = 



output sequence v assuming (74) 
no signal is present ^ 

for Gaussian-stationary detector noise. Since the detector is linear, the probability 

/ probability of observing 
P{y\}i) = output sequence v assuming (75) 
signal h is present y 

is just 

Piv\h)^Piv-Vh\0), (76) 

where is the detector response to the gravitational wave signal h. The ratio of 
these two probabilities, 

Mv\h) ^ (77) 



p{v\oy 

termed the likelihood function, is the odds that the data v is a combination of 
signal Vh and noise, as opposed to a noise alone. For a given observation v the 
likelihood can be viewed as a function of hypothesized signal h, in which case it 
has a convenient interpretation in terms of plausibility: in particular, A{v\h) can 
be interpreted as the plausibility that the signal h is present given the particular 
observation v. (The hkelihood is not, however, a probabihty.) This meaning of the 
hkehhood is independent of the statistical character of the noise. The difficulty, 
if the noise is not Gaussian-stationary, is in evaluating A. 



3.5.5 The two-time correlation function 

The correlation function C2(r) describes the statistical relationship between pairs 
of samples drawn from the random process n{t) at times separated by an interval 
r. Given two samples separated in time by r, a non-zero correlation C2(t) cor- 
responds to an increased ability to predict the value of one member of the pair 
given the other. 

The correlation function C2(r) is bounded by ±6*2(0), suggesting that we 
define the correlation coefficient 



R,{t) = C2{t)/C2{0), 



(78) 



which is bounded by ±1. If the correlation coefficient is zero for some r, then 
samples taken an interval r apart are entirely uncorrelated: knowledge of one 
does not lead to any increased ability to predict the other. A positive correla- 
tion coefficient tells us that the two samples are more likely close to each other 
in magnitude and sign than not, while a negative correlation coefficient tells us 
that the two samples are likely close to each other in magnitude but of opposite 
sign. The larger the coefficient magnitude the greater the tendency. When the 
correlation coefficient is unity then the correlation is perfect: i.e., when it is +1 
the two samples are always equal, and when it is —1 the two samples are always 
of equal magnitude but opposite sign. 

3.5.6 Noise Power Spectral Density 

Consider for a moment a simple harmonic oscillator — e.g., a pendulum — coupled 
weakly to a heat bath. The heat bath excites the oscillator so that its mean energy 
is ksT. Since the coupling to the heat bath is weak, the phase of the oscillator 
progresses nearly uniformly in time with rate cuq corresponding to the oscillator's 
natural angular frequency. Over long periods, however, the continual, random 
excitations of the oscillator cause the phase to drift in a random manner from 
constant rate. 

Now suppose that we sample the position coordinate of the oscillator at inter- 
vals separated by exactly one period 2t: /ujq. Since the coupling to the heat bath 
is weak the samples are very nearly identical: in fact, were it not for the contact 
with the heat bath, they would be exactly identical. Thus, we expect that the 
correlation coefficient corresponding to an interval equal to an oscillator period 
should be nearly unity. Continuing to focus on samples taken at intervals equal 
to exact multiples of the period, we expect that the correlation coefficients should 
remain large for small multiples, but should decrease as the interval increases 
since contact with the heat bath will lead, as time increases, to greater drift in 
the phase. 

On the other hand, suppose that we sample the position coordinate of the 
oscillator at intervals separated by exactly odd integer multiples of a /laZfperiod 
"k/ujq. Now we expect the correlation coefficient to be nearly equal to —1 for small 
intervals, decreasing in magnitude to as the interval increases. 

Contact with a heat bath can take place in many ways, leading to subtly differ- 



ent correlation functions. Figure ^ shows the correlation function corresponding 
to two different kinds of heat bath contact: that which leads to velocity damping 
and that which leads to structural damping.!^ Note how, pictured in this way, 
there is apparently little difference between these two damping measures. 

Since the correlation function is so oscillatory we are immediately led to con- 
sider its Fourier transform. In this case, since C2(t) is an even function of the lag 
r, we consider the cosine transform, which we term the one-sided power spectral 
density: 



{One-sided refers to the fact that, in choosing a cosine transform, we have ef- 
fectively folded the power in negative frequencies into the power at positive fre- 
quencies; so, the Sv{f) includes the power at frequencies whose magnitude is 
|/|.) Figure ^ shows the power spectral densities corresponding to the correlation 
functions of figure |^. The strongly oscillatory nature of these functions shows 
up as a large peak at the oscillator resonant frequency (normalized to unity). In 
addition, however, the PSD shows clearly the very different off-resonance char- 
acter of the noise. Noise from a structurally damped system rises in amplitude 
as the frequency falls below resonance, unlike the noise contribution from a vis- 
cously damped system; similarly, noise from a structurally damped system falls 
more steeply with frequency above resonance than does the noise from a viscously 
damped system. To see the same in the correlation function would require close 
inspection of the trends of the correlation function envelopes over very long lags. 

Thus, even though it is completely equivalent to the correlation function, the 
power spectral density is often a more useful characterization of the noise char- 
acter. In the case of Gaussian noise its equivalence to the correlation function 
guarantees it is also a full characterization of the detector noise. When the noise 
is not Gaussian, there are analogous spectra associated with the higher order cor- 
relations: for example, the bispectrum is the 2-dimensional Fourier transform of 



and so on. These higher order spectra and their magnitudes play the same role 
for the higher-order correlation functions as the power spectral density plays for 
the auto-correlation function. 




(79) 
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Fig. 3. The correlation functions for a harmonic oscillator in contact with a heat 
bath. Contact with a heat bath leads to damping; the nature and degree of the 
contact determine the character of the system's noise. In this figure we show 
the correlation function, over several periods, for two different kinds of heat bath 
contact with the same on-resonance damping. The upper panel corresponds to 
a viscous damped harmonic oscillator; the lower panel corresponds to a struc- 
turally damped oscillator. The difference between the two correlation functions is 
apparently very subtle. 
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Fig. 4. The power spectral density of the two processes whose correlation functions 
are shown in figure |^. Note that, while the correlation functions appear very 
similar as functions of time, strong differences show up in the power spectral 
densities as functions of frequency. 



3.6 Signal-to-noise ratio 

When is a gravitational wave "detectable" ? We haven't yet explored the meaning 
of "detection" qualitatively, let alone quantitatively; nevertheless, we have an 
intuitive feeling that a signal ought to be detectable if the detector's response to 
the signal is greater than the intrinsic noise amplitude. Let's develop that idea a 
bit. 

Suppose that we have a detector with noise power spectral density <S'„(/) and 
particular output v{t), which consists of a signal Vh{t) superposed with detector 
noise Vn{t). The variance of v{t), over an interval [0, T], is 

= y^^dtvitr (81) 

The noise is a random process; so, then, is a^. Focus on the ensemble average of 
and look in the frequency domain: 

f ^ pilW (83) 

2 



{\Vn{f)f+\Vk{f)\") (84) 



T 

where the final equality follows when we recognize that the noise is independent of 
the signal. The contribution to the mean signal variance thus consists of separate 
contributions from the signal and from the noise. 
The ratio 



\Vn{f)f 



evidently tells us which — signal or noise — is expected (note ensemble average!) 
to contribute more to the amplitude of the detector output in a unit bandwidth 
about frequency /. We can compute a similar, dimensionless quantity over the 
full bandwidth: 

,f^m^=,r,fM]i (86) 

• \Vn{f)\/T Jo S^if) 

tells us which of the signal Vh or the noise v„ is expected to contribute more to 
the variance of the output v. 

Given a particular sample of detector output v we don't know, a priori, what 
part is Vn and what part (if any) is Vh. Consider a quantity that we can calculate 



J — ( 



directly from the detector output v: 



?2 = 4 



df 
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The integrand is evidently the ratio of the actual contribution to the signal vari- 
ance in a unit band about frequency / to the contribution that would be expected, 
in the same band, from noise alone. Not surprisingly, the ensemble mean jg 



1 + 4 



df 



We refer to as the signal-to-noise ratio, or SNR.[tt] 

Our construction of p^ has been physically motivated. It turns out, however, 
that exactly this same quantity arises from a consideration of the probability 



P(f |0), which we explored in § |3.5.3| . In that section we found, for Gaussian- 
stationary noise, 

'-k^f,k2ov[j]\\C-%,v[k] 
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With just a little algebra, however, the argument of the exponential can be rewrit- 
ten as0 
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(90) 



where the periodic sequence g[j] is related to the discrete Fourier transform of the 
sequence g[k]: 

N-l 

9[j] = E 

k=-{N-l) 

and V[k] is just v[k] zero-padded for negative k: 
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(91) 
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These summations can be regarded as approximations to integrals, in which case 
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II Note that this definition of is different, by the additive factor of unity, than used elsewhere 
in the gravitational wave literature. 




Frequency 

Fig. 5. An imagined gravitational wave signal (upper-left panel), detector noise 
(middle-left panel) and total detector output (signal+noise, lower-left panel). 
Note how the signal is not evident to the eye in the detector output. The right- 
hand panels show the power spectra of the corresponding left-hand panels; again, 
the signal power is not evident relative to the noise power. 

Hence, the SNR associated with the observed detector output v is closely related 
to the probability that f is a sample of just detector noise, with no gravitational- 
wave signal present. The larger p^, the smaller this probability. Should we observe 
detector output with large p^, then, we are not too far wrong to be suspicious that 
we have seen evidence for gravitational waves. 

To make this last judgment — which involves making quantitative the notion 
of "large" p^ — we need to know the probability distributions of the SNR in both 
the presence and absence of a signal: after all, since noise is a random process 
there is some non-zero probability that, in any given observation, p^ will take on 
any particular value, large or small. We return to consider this point in §^ 



3.6.1 Matched Filtering 



Calculating defined by equation |8^ does not require or make use of any infor- 
mation about the gravitational radiation source. Suppose that we know, a priori, 
the radiation waveform has the shape Vh{t), and that the question is whether the 
corresponding signal aVhit — to), for some unknown constants a and to, is present 
in the detector observed output v (t) . Can we make use of this information — the 
signal shape Vh{t) — to boost our ability to observe the signal? 

The answer is yes. To illustrate, figure || shows an imagined Vh, f « and v equal 
to Vh + Vn. in the left-hand panels, and the corresponding power spectra in the 
right-hand panels. For this illustration we have assumed that the noise is white 
across the detector bandwidth. The signal is not apparent to the eye in either v 
or its power spectrum Pv{f)- Figure |] shows, in the top panel, the filter output 
when just Vh is passed through the filter K with impulse- response Vh set equal to 

Vh- 



where 



v'(t) = / dTv{T)K{t-r) (94) 

= vUt)+v',{t) (95) 

K{t) = Vh{t) (96) 
t 



v'^{t) = / dTVn{r)Vh{t + r) (97) 

J— oo 

V'hit) = f dTVh{T)v'^{t + T). (98) 



Without loss of generality we assume is non-zero only for positive t. The 
filtered detector output v'{t) consists of a signal contribution f^(t) and a noise 
contribution f^(t). These are shown in the top and middle panels of figure H, 
respectively. The bottom panel of figure || shows the filter output v' (equal to 
v'h + ^n)- The presence of the "signal" is now much more evident. 

The filter we have chosen has reduced the total power in the noise relative to 
that in the signal. How it does this is apparent by considering the power spectra 
in figures ^ and ^. In figure ^, the power in Vh is seen to be confined to a very 
narrow bandwidth about the frequency of the damped sinusoid. At its peak the 
signal power is about 5 dB greater than noise power. Nevertheless, the total noise 
power, integrated over the full bandwidth, is much greater than the signal power 




Fig. 6. The output of the filter described in equation 95 when just the signal Vh 
is filtered (upper panel) and when the detector output, consisting of signal and 
noise, is filtered (lower panel). In contrast to the lower-right panel of figure ^, the 
"signal" {i.e., the upper panel) is quite evident even in the presence of noise. 



and, consequently, the signal is overwhelmed by the noise (cf. the bottom panel 
of figure 

Now consider v' . The filter applied to the signal has the impulse response 
of the signal, or the squared magnitude frequency response given by the power 
spectrum in the top panel of figure ^. This is matched to the signal, in the sense 
that the power passed is in the band where the signal power is large and the power 
stopped is in the band where the signal power is small. Thus, what survives in 
v' is the signal power, together with only that noise power in the narrow band 
where the signal power is large. The signal to noise of the filtered detector output 
v' is correspondingly much higher in the presence of the signal than is the signal 
to noise ratio of v. 

This example is illustrative. In fact, we can ask, for an arbitrary signal Vh em- 
bedded in noise with power spectrum S^i^f), for the linear filter that maximizes 



the ratio of the mean-square signal contribution to the mean-square noise contri- 
bution. That filter is referred to as the Wiener matched filter; in the frequency 
domain and for weak signals it is (up to an overall constant) 

m - f§r. (99) 

More generally, additional information is always useful for increasing the our 
ability to detect a signal. This is true even that information is not as complete as 
knowing the waveform. For example, consider the case where we know the signal 
spectrum, but not its waveform. In the frequency domain, we thus know the 
signal amplitude at each frequency, but not the corresponding phase. In the case 
where the waveform is known, we constructed the filter making full knowledge 
of both amplitude and phase information. We can also construct a filter that 
passes power in a given bandwidth, without regard to its phase. This filter will 
emphasize power in the bands where the ratio of signal power to noise power is 
relatively large over bands where the ratio is small; consequently, it will increase 
our ability to detect a signal whose spectrum is known in the same way that a 
matched filter increases our ability to detect a signal whose waveform is known. 

3.7 The effective noise power spectral density 

How does one compare different detectors, with different response functions and 
different noise power spectral densities? 

One possibility is the "performance benchmark" : choose a prototypical source, 
evaluate the signal-to-noise that the source would give in the different detectors, 
and determine finally which detector is most likely to observe the source at a given 
level of confidence.ii 

This kind of judgment depends critically on the source: using different sources 
as your benchmark can lead to different conclusions. For example, sources whose 
power is concentrated at different frequencies focus attention on the detector 
noise at those frequencies. Thus, while benchmarking detectors against partic- 
ular sources can be a powerful tool for comparing their relative performance, it 
is also a tool with a very narrow focus. We need some other way to compare the 
capability of detectors with a less specific emphasis on source. 

An important tool for making this more general comparison is the effective 
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Fig. 7. The power spectral density of an effective stochastic gravitational wave 
signal that would mimic the noise in the output of a modern bar detector. Plotted 



is yShif) vs. frequency /. 
power spectral density Sh{f), 



SM) ^ (100) 



where R{f) is the detector response function and is the detector noise power 

spectral density. The quantity Sh{f) describes an effective detector noise: it is the 
power spectral density of a stochastic gravitational wave signal that would have to 
be applied to a noise-free detector in order that the corresponding response have 
power spectral density S^if). Over frequency bands where Sh{f) is small, the 
detector is relatively sensitive; over frequency bands where it is large, the detector 
is relatively insensitive. 

The effective power spectral density Sh{f) has both a source and detector in- 
dependent meaning, making it a particularly useful quantity for comparing gravi- 
tational wave detectors or for comparing a detector to a source. With it, one can 



rank detectors according to their overall noise in a given bandwidth, e.g., 

F.(/i,/2)= f^'dfs.if), (101) 



/l 

or define an effective band (/o — A//2, /o + A//2) over which the detector has 
greatest sensitivity, e.g., 

r _ Io^dff/S,{f) 

= iTdf/s,U) ^ ^ ^ 

Finally, since the noise is referred directly to the amplitude of incident gravi- 
tational radiation, one can calculate the expected SNR of a given signal in the 
detector without reference to the detector's response function: 

~ 2 

roo 

p2 = 1+4/ df 

Jo 



R{f)h{f) 



s.if) 

~ 2 
roo h(f) 

= 1 + 4/ rf/^TTT- (104) 

Figure |^ shows the modeled Sh{f) for a modern bar detector, while figure 
^ shows Sh{f) for a model of the first-generation LIGO instrumentation. Note 
how the bar detector noise is particularly small in two narrow bandsQ about 
the resonant frequencies of the two mode system consisting of the bar and its 
transducer, while the interferometer achieves its peak sensitivity over a much 
broader bandwidth. 

3.7.1 An aside: noise in bar detectors 

It is a common misconception that bar detectors are intrinsically narrow-band 
detectors. While the amplitude of a resonant detector's response is greatest for 
signal power in the neighborhood of the resonance, the thermal excitation of the 
bar is also concentrated in this band as well. The net result is that the contribution 
of the bar's thermal noise to the power spectral density expressed in units of /i^/Hz 
is effectively independent of frequency. 



''Since the bar detectors "sensitivity" 1/Sh is niulti- modal it is more appropriate to define the 



effective band, as in eq. f02 and 103, separately about each peak. 



Fig. 8. The power spectral density of an effective stochastic gravitational wave 
signal that would mimic the noise in the output of first generation LIGO instru- 



To understand how resonant detectors become narrow band instruments, con- 
sider how the signal appears in the electronics that follow the transducer. The 
resonant character of the detector leads to large amplitude motion for signal power 
near the resonant frequency and small amplitude motion for signal power far from 
the resonance. Correspondingly, the amplified signal is large near to, and small 
far from, the resonance. The amplifier contributes its own noise, however, which 
is approximately white at the amplifier output. Thus, compared to the signal pre- 
sented for amplification, the amplifier noise is relatively large far from resonance 
and relatively small near to resonance. 

In present-day resonant cryogenic detectors the bandwidth is limited by am- 
plifier noise, referred back to h through the response function. 

Since it is the amplifier noise, when referred back to h through the response 
function of the resonant bar, that limits the instrument bandwidth, why make the 
bar resonant at all? The purpose of making the detector resonant is to provide 



mentation. Plotted is J Sh{f) vs. frequency /. 




mechanical amplification for the signal, so that, at least in a narrow bandwidth, 
it is much stronger then the limiting noise source: amplifier noise. Signal power 
at or near resonance leads to a large excitation of the bar, which translates into a 
large input to the amplifier; thus, the resonance of the bar amplifies the incident 
gravitational wave signal relative to all the noise sources that follow, including the 
limiting amplifier noise. 

3.8 Conclusion 

Gravitational wave detectors are characterized by their antenna patterns, which 
describe their differential sensitivity to radiation incident from different directions 
and with different polarizations, their response functions, which describe the dif- 
ferential amplitude of their response to signals of different frequency, and the 
character of their noise. 

The distinction between the response function and the antenna pattern is 
sometimes an artificial one: the response function can (and for interferometric 
detectors, does) depend on the incident direction of the radiation. 

Much of the experimental craft is devoted to making the detector noise ap- 
proximately stationary and Gaussian (or in making the signal so large that the 
character of the noise is not significant for measurements of interesting precision) . 
Stationary noise can be characterized by evaluating the correlations among sam- 
ples taken at different relative times. For Gaussian noise, all the correlations are 
known once the pairwise correlation function is measured. 

While the correlation functions are good conceptual tools for understanding 
the character of stationary detector noise, a more useful tool, fully equivalent, 
is the noise spectrum: the Fourier transform of the correlation function on its 
time arguments. The noise power spectral density, in particular, is a particu- 
larly important and useful tool for characterizing the disposition of detector noise 
power. 

4 Characterizing detection 

What does "detection" mean? Let's try to frame an answer by posing a specific 
question — e.g., "with what confidence can we conclude that, in the last hour, 
gravitational waves from a new core collapse supernova in the Virgo cluster of 



galaxies passed through our gravitational wave detector?" — and exploring its 
meaning. 

It turns out that, straightforward as it seems, there arc two different ways of 
interpreting this question; correspondingly, there are two different meanings that 
"detection" can take. How we mean the question — or the k;ind of answer that we 
want to tak;c away — determines the kind of analysis that we need to undertake 
with data collected at a gravitational wave detector. 

4.1 Learning From Observation 

"With what confidence can we conclude that, in the last hour, gravitational waves 
from a new core collapse supernova in the Virgo cluster of galaxies passed through 
our gravitational wave detector?" 

As with most questions of detection, even before examining our observations we 
have some expectation of the answer. In this case, we know the rate of supernovae 
and this leads us to expect, on average, one such core collapse every 4 months; 
consequently, we believe the probability is approximately 3.4 x 10~^ that in any 
given hour — including the last — gravitational waves from a new Virgo cluster 
supernova were incident on our detector. 

Probability, as we have used it here, means degree of belief. In this instance, 
our degree of belief coincided with the expected frequency of supernova events. 
This need not always be the case: we can assess degree of belief even when we 
can't assess relative frequency. For example, suppose that I have a coin that is 
known to be heavily biased toward either heads or tails. What is your degree 
of behef that, when I next flip the coin, it will land heads-up? Without telling 
you the direction or amount of the bias, you can't evaluate the expected relative 
frequency of heads or tails. You can, however, quantify your degree of belief: 
having no more reason to believe that the bias is toward heads than towards 
tails, you have no more reason to believe that the coin will, when next flipped, 
land heads-up than that it will land heads-down. Your degree of belief in either 
alternative, then, is 1/2. 

One does not have to search either long or hard to find examples from, e.g., 
astrophysics, where probability as "degree of behef" exists and probabihty as 
"expected frequency" does not. For example, what is the probability that there 
exists a cosmological stochastic gravitational wave signal with a given amplitude 



and spectrum? In this case, "expected frequency" has no meaning: there is only 
one Universe, and it either does or does not have a stochastic gravitational wave 
background of given spectrum and amplitude. 

After we examine the output of our gravitational-wave detector, our degree 
of belief in the supernova proposition may change: we may, on the basis of the 
observations, become more or less certain that radiation from a supernova passed 
through our detector. How do observations change our degree of belief in the 
different alternatives? 

To explore how our degree of belief evolves with the examination of observa- 
tions we need to introduce some notation: 

^ proposition that gravitational waves from a 
Hq — new supernova in the Virgo cluster did not , (105) 
pass through our detector in the last hour j 

(our prior knowledge of astrophysics, including | (^io6) 
our best assessment of the supernova rate I 

g = observations from our gravitational wave detector ) , (107) 

P{A\B) = degree of belief in A assuming that B is true ^ , (108) 

-I A — (logical negation of proposition ^4) . (109) 

In this notation, P{Ho\I) is the degree of behef we ascribe to the proposition 
that no gravitational waves from a core collapse supernova in the Virgo cluster 
passed through our detector in the last hour, given only our prior understanding 
of astrophysics; similarly, P{Ho\g,X) is the degree of belief we ascribe to the 
same proposition, give both the observation g and our prior understanding of 
astrophysics. 

To understand how P{Ho\I) and P{Ho\g,I) are related to each other we need 
to recall two properties of probability. The first is unitarity: probability summed 
over all alternatives is equal to one. In our example, the two alternatives are that 
a supernova occurred or it did not: 

P{Ho\g,I) + PhHo\g,I)^l- (HO) 

The second property we need to recall is Bayes Law, which describes how condi- 
tional probabihties "factor" : 



P{A\B, C)P{B\C) = P{A, B\C) = P{B\A, C)P{A\C). (Ill) 



Combining unitarity and Bayes Law it is straightforward to show that 
where 

A{g) = P{g\^Ho,I)/P{g\Ho,I) (113) 
, , / probabihty that g is a sample of \ , , 

Pig\Ho,l) = \ y y (114) 

\ detector output when Hq is true / 

P{g\^Ho I) = I P^°^^^^^^*y 5- is a sample of \ ^^^^^ 
\ detector output when Hq is false / 

The two probabilities P{g\HQ,X) and P{g\^Ho,X) depend on the statistical 
properties of the detector noise and the detector response to the gravitational wave 
signal. In some cases they can be calculated analytically; in other circumstances 
it may be necessary to evaluate them using, e.g., Monte Carlo numerical methods. 
Regardless of how one approaches data analysis the detector must be sufficiently 
well characterized that these or equivalent quantities are calculable. 



Equation |112| describes how our degree of belief in the proposition -^Hq evolves 
as we review the observations. If A is large compared to the ratio P{Hq\X) to 
P{-iHq\I) then our confidence in -li^o increases; alternatively, if it is small, then 
our confidence in -^Hq decreases. If A is equal to unity — i.e., the observation g 
is equally likely given Hq or -lifg — then the posterior probability P{Ho\g,I) is 
equal to the prior probability P{Hq\I) and our degree of belief in Hq is unchanged: 
we learn nothing from the observation. 

We can now answer the question that began this section. We understand 
confidence to mean degree of belief in the proposition that radiation originating 
from a new supernova in the Virgo cluster was incident on a particular detector 
during a particular hour. In response we make a quantitative assessment of our 
degree of belief in that proposition — the probability that the proposition is true. 



4.2 Guessing Natures State 

Begin again: 'With what confidence can we conclude that, in the last hour, 
gravitational waves from a new core collapse supernova in the Virgo cluster of 
galaxies passed through our gravitational wave detector?" 



As before, we have the hypothesis Hq and its logical negation, -iHo. The 
gravitational waves from a new Virgo cluster supernova either passed through our 
detector, or they did not. Our goal is to determine, as best we can, which of these 
two alternatives correctly describes what happened. 

We decide which alternative is correct by consulting our observation g. Op- 
erationally, we adopt a rule or a procedure that, when applied to g, leads us to 
accept or reject Hq. The question that began this section asks us for our degree 
of confidence in the most reliable rule or procedure. 

There are many procedures that we can choose from. Some are just plain 
silly: for example, always rejecting Hq is a procedure. Similarly, accepting Hq if 
a flipped coin lands heads is a procedure. Some procedures are more sensible: we 
can calculate a characteristic amplitude from the observation {e.g., a signal-to- 
noise ratio) and reject Hq if the amplitude exceeds a threshold. Nature doesn't 
always speak clearly; additionally, some crucial information is often hidden from 
us. Consequently, no procedure will, in the end, be perfect and every rule will, on 
unpredictable occasions, lead us to erroneous conclusions. Still, some procedures 
are clearly better than others: the question is, how do we distinguish between 
them quantitatively? 

Better procedures are those that are less prone to error. Consequently, we focus 
on the frequency with which different procedures err. For our simple problem, 
where we want to decide only if we have or have not observed the radiation from a 
supernova (reject or accept Hq), there are two kinds of errors a decision procedure 
can make: 

1. If no radiation is present {Hq true), the rule may incorrectly lead us to 
conclude that radiation is present: a false alarm, or type I, error. 

2. If radiation is present {Hq false), the rule may incorrectly lead us to conclude 
that radiation is absent: a false dismissal, or type II, error. 

The false alarm frequency is generally denoted a while the false dismissal fre- 
quency is denoted /3. 

If we have an ensemble of identical detectors, each observing simultaneously 
the same system for which Hq (or -iHq) is true, and we apply our rule to each 
observation, then the fraction of observations in the ensemble that lead us to 
reject (accept) Hq is just the false alarm (dismissal) frequency. False alarm and 



false dismissal frequencies can be interpreted as probabilities: in particular, the 
probability of our rule giving an incorrect result. 

Even in the simple case at hand (a single hypothesis that we must accept or 
reject), there are at least two distinct kinds of errors that an inference procedure 
can make. Our measure of a rule's reliability thus involves at least two dimensions, 
and may involve more. How, then, do we order rules to settle upon a best, or 
optimal, rule? 

To rank rules we must reduce the several error measures that describe a proce- 
dure's performance to a single figure of merit. How we choose to do this depends 
on the nature of our problem. In our case, rules that distinguish between Hq and 
-li^o are characterized by their false alarm and false dismissal frequencies; conse- 
quently, our criteria for ranking rules should depend on our relative intolerance to 
false alarms and false dismissals. For example, if we are testing for the presence 
of antibodies in an effort to diagnose and treat a serious illness, we might be very 
concerned to keep the false dismissal rate low, and not nearly as worried about a 
high false alarm rate: after all, a false dismissal might result in death, while a false 
alarm only in an unnecessary treatment with less serious repercussions. Judges or 
juries in criminal trials faces different concerns: false dismissals let criminals go 
free, while false alarms send the innocent to prison — neither alternative being 
very palatable. Finally, in the case of gravitational wave detection, we may (at 
least initially) be very concerned to avoid false alarms, even at the risk of falsely 
dismissing many real signals. 

Thus, in order to provide a relative ranking of different inference procedures 
for detection or parameter estimation we must construct an ad hoc figure of merit 
that refiects our sensitivity to an incorrect decision. We term the best rule, under 
that ad hoc criteria, the "optimal" rule. "Optimality" , however, is a relative 
concept: if the criteria change, the "optimal" rule changes also. In the three 
examples given above, the criteria might be 

• medical diagnosis: fix a maximum acceptable false dismissal rate and choose 
the rule that, among all rules whose false dismissal rate is so constrained, 
has the minimum false alarm rate; 

• criminal justice: choose a rule whose weighted total error a cos (f) + /S sin (j) 
is minimized (0 being a matter of personal choice for an individual judge or 
juror); 



• gravitational wave detection: fix a maximum acceptable false alarm rate and 
choose the rule that, among all rules whose false alarm rate is so constrained, 
has the minimum false dismissal rate. 

False alarm and dismissal rates describe our confidence in the long-run behav- 
ior of the associated decision rule. To understand the implications of this measure 
of confidence, suppose that we have not one, but N independent and identical 
detectors all observing during the same hour. We use the same test, with false 
alarm rate a and false dismissal rate (3, on the observations made at each detector, 
and find that, of these N observations, m lead us (through our inference rule) to 
reject Hq and N — m lead us to accept Hq. For a concrete example, suppose a is 
1%, N is ten and m is three. 

The probability of obtaining this outcome when the signal is absent {Hq is 
true) is the probability of obtaining m false alarms in N trials, or 

P(m|//o, N) = - af-"^. (116) 

In our example, P(m|ifo) ^) evaluates to 1.1 x 10"*^. It is thus very unlikely that 
we would have made this observation if the signal were absent. Does this mean 
we should conclude the signal is present with, say, 99.99% confidence? 

No! P{m\Ho, N) describes the probability of observing m false alarms out 
of N observations. When the signal is present, however {i.e., when Hq is false), 
there are no false alarms and both a and P{m\HQ, N) are irrelevant. There are, 
however, N — m false dismissals; thus, the relevant quantity is P{m\-iHo, N), the 
probability of observing N — m false dismissals: 

If, in our example, the false dismissal rate j3 is 10%, then the probability of 
observing seven false dismissals out of ten trials is is 8.7 x 10~^. 

The particular outcome of our example — three positive results out of ten trials 
— is, in the grand scheme of things, very unlikely; nevertheless, what is important 
to us is that it is more unlikely to have occurred when the signal is present then 
when it is absent. Despite the apparently overwhelming improbability of three 
false alarms in ten trials, it is nevertheless, slightly more likely than the alternative 
of seven false dismissals in ten trials. 



We can now answer the question that began this section. We understand 
that question to ask for the error rate of the best general procedure for deciding 
between the alternative hypotheses. There is an implicit assumption regarding 
the decision criteria, which tells us what "best" means in this context. In the 
context of these criteria, we calculate the error rates for different inference rules, 
choose rank the different rules, and find the best rule and its corresponding error 
rates. 

Contrast this with our understanding of the identically worded question as 
we understood it in the previous subsection. There, we understood confidence to 
mean the degree of belief that we should ascribe to alternative hypotheses; here, we 
understand confidence to refer to the overall reliability of our inference procedure. 
There we responded with a quantitative assessment of our degree of belief in 
the alternative hypotheses, given a particular observation made in a particular 
detector over a particular period of time; here we responded with an assessment 
of the relative frequency with which our rule errs given each alternative hypothesis. 
There we did not make a choice between alternative hypotheses; rather, we rated 
them as more or less likely to be true in the face of a particular observation. Here, 
on the other hand, we do make choices and our concern is with the error rate of 
our procedure for choosing, averaged over many different observations and many 
different decisions. 

Analyses like the ones in this section, where probability is interpreted as the 
limiting frequency of repeatable events and the focus is on false alarm and false 
dismissal frequencies, are termed Frequentist analyses. They have particular util- 
ity when it is possible to make repeated observations on identical systems: e.g., 
particle collisions in an accelerator, where each interaction of particle bunches 
is a separate "experiment." Analyses like those in the previous section, where 
probability is interpreted as degree-of-belief and the focus is on the probability 
of different hypotheses conditioned on the observed data, are termed Bayesian 
analyses. Bayesian analyses are particularly appropriate when the observations 
or experiments are non-repeatable: e.g., when the sources are, like supernovae, 
non-identical and destroy themselves in the process of creating the signal. In this 
case we are interested in the properties of the individual systems and would prefer 
a measure of the relative degree of belief that we should ascribe to, for example, 
the proposition that the signal originated from a particular point in the sky. 

That Bayesian and Frequentist analyses are different does not imply that one is 



right and the other wrong. Bayesian and Frequentist analyses do not address the 
same questions; so, they are not required to reach "identical" conclusions. On the 
other hand, it may well be that one analysis is more appropriate or responsive to 
our concerns than the other. We can only make the choice of appropriate analysis 
tools when we understand the distinction between them. 

5 Gravitational Radiation Sources 

In this section we review briefly some of the different kinds of sources that are, at 
this writing, thought to be "important" for the generation of large interferometric 
detectors now under construction. "Important" is a term that requires definition 
in this context. Clearly, sources that we don't expect to detect, or to be able to 
detect, are unimportant. Detectable sources must radiate significant energy in 
the bandwidth where these detectors are most sensitive (ranging from the tens to 
hundreds of Hz) over a reasonable observation period. For periodic sources, this 
is the integrated power over a period of several months to perhaps as much as one 
to two years; for burst sources, this means that the expected rate of detectable 
bursts must be at least several per year. 

The assessment of source strength and number or rate is difficult to make for 
most sources. Very often the radiation strength depends on physics and astro- 
physics that we don't know or understand in the requisite detail. For all proposed 
sources the rate, number or even existence of sources whose signal strengths are 
large enough to be detectable is difficult to ascertain. This is not surprising: what 
we know of the heavens we know principally through electromagnetic observa- 
tions; however, it is in the nature of gravitational wave sources that they leave 
little electromagnetic evidence of their existence. 

Finally, since we are, with these instruments, looking at the universe in a 
fundamentally new way, we must keep an open mind to the possibility of sources 
unimagined: in this I side with John Haldane, who said (in a different context) 
"My own suspicion is that the universe is not only stranger than we suppose, but 
stranger than we can suppose." 

One final note: until now, we have been careful to keep all factors G and c 
in our expressions for gravitational effects. Here and henceforth, we will write all 
expressions in units where G and c are unity: e.g., units of length, with conver- 
sion factors from grams and ergs to centimeters as given in equations |^ and ^. 



These conversion factors can be invoked to find expressions in terms of quantities 
expressed in more conventional units. 



5.1 Burst sources 

5.1.1 Compact binary inspiral 

The source most-talked-about for the interferometric detectors now under con- 
struction are binary systems consisting of two compact, stellar mass objects — 
either neutron stars or stellar mass black holes. Like a rotating dumbbell, a binary 
star system has a large, accelerating quadrupole moment, which makes it (for its 
mass) a strong gravitational radiation source. The radiation carries away orbital 
binding energy and orbital angular momentum, which leads to a faster and more 
compact orbit. Kepler's Third Law relates the orbital frequency /orb' semi-major 
axis a and total system mass M by 

/orb = (118) 
Consequently, the radiated power is, in order of magnitude. 
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Note that radiated power increases as the orbit decays: 
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For binaries that can become sufficiently compact the power radiated gravitation- 
ally will, in the end, become large enough to dominate the system's evolution. 
Since the radiated power increases as the orbit decays, the system will then decay 
at an ever increasing rate, with ever increasing radiation amplitude, frequency 
and power, until the components coalesce. It is the radiation from this inspiral, 
for binary systems of neutron stars or black holes, that is seen as an important 
source for the LIGO and VIRGO detectors. 

Why compact components, like neutron stars or black holes, and why stellar 
mass, and not more or less massive? Recall that the proposed interferometric de- 
tectors have their greatest sensitivity at approximately 150 Hz. The quadrupole 
radiation from a binary system is at twice the system's orbital frequency; corre- 
spondingly, if the radiation is to be in the bandwidth of these detectors, the binary 



systems themselves must exist with orbital frequencies of at least 75 Hz. Kepler's 
Third Law places a lower bound for us on the matter density of the components, 
which must be much greater than the total system mass divided by the cube of 
the orbital radius: 

M 

Pi,2 » — (121) 
^ (^r/gw)' (122) 
~ 1.5 X 10^2 ( Jm_] gW. (123) 

V 100 Hz; ^' ^ ' 

Thus, irrespective of the total system mass, if a binary system is to radiate in a 
band where these detectors are sensitive the central density of its components can- 
not be much less than nuclear density. With this we are forced, for astrophysical 
objects, to restrict attention to neutron stars or black holes. 

The nuclear and super-nuclear equation of state place upper limits on the 
neutron star mass, which does not apply for a black hole. The dynamics of 
the binary orbit, however, does place an upper limit on the mass of the black 
hole binaries that the ground-based interferometric detectors may observe. With 
every orbit the binary radiates away more of its binding energy, leading to a more 
compact orbit. Eventually the system coalesces: the two components merge, 
collide, or tidally disrupt. Even if we imagine that the components are point 
masses, so that there is no tidal disruption or collision that would terminate 
the inspiral signal at some finite orbital frequency, relativity appears to impose 
a maximum orbital frequency on binary systems. For approximately symmetric 
binary star systems {i.e., those with equal mass components) this limit iJ^l 

/max ^710^:^ Hz, (124) 

where M is the system's total mass. Thus, the component black hole masses must 
be less than 15 if the inspiral signal is to survive into the bandwidth where 
the detector is most sensitive. 

It is currently thought that, during the epoch when the radiation from the 
binary is in the bandwidth where the LIGO and VIRGO detector sensitivity is 
greatest, the binary components are well approximated as point masses for the 
purpose of computing the radiation and orbital evolutiorS (There is some small 
suggestion that resonant tidal interactions may complicate this picture.il) Dur- 
ing this epoch, the gravitational fields that determine the binary evolution are 



sufficiently strong that ffist order perturbation theory is not adequate to compute 
the orbits; nevertheless, the fields are not so strong that computing the orbits 
and the radiation via higher order perturbation theory is impractical.00 For 
this overview, no additional insight is gained by considering anything higher the 
quadrupole formula radiation, in which case the excitation of the detector — an 
effective h{t) that is a superposition of the radiation in the two polarization states 
of the wave — isSi^ 

h{t) = {nfMf^ cos (125) 



where 
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= l'2nf(t')dt', (129) 

is the cosmological luminosity distance to the source, mi and m2 are the binary 
system's component masses, z is the source's cosmological redshift, l is the angle 
between the binary's angular momentum axis and the line of sight to the detector, 
and To is a constant of integration. 

What can we determine through observation of the signal from such a system? 
The signal-to-noise, of course, which takes on a particularly simple formU'El: 

where tq is a characteristic distance that depends only on the effective power 
spectral density of the source, 

fGM^\' 5 f 243 r df {ttGM^/c')' 



(127) 
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and the average denoted by the over-bar is over both an ensemble of detectors 
and all relative orientations of the source and the detector. For the initial LIGO 
and VIRGO detectors, is about 13 Mpc. In order that we are confident that we 
have seen a source, the SNR should not be much less than about 65 in a single 



detector; so, we don't expect to see sources from distances beyond more than a 
few tq. 

How many of these sources can LIGO expect to see? Unfortunately, we know 
very httle about the rate of compact binary coalescence, except that it is rare. 
Black-hole/black- hole binaries are, of course, invisible to us except through gravi- 
tational waves. Binaries involving neutron star component (s) are observable to us 
only if one of the components is a pulsar. Pulsars are observable only if they are 
not too distant (in our galaxy or its satellite globular clusters) and if the pulsar 
beam intersects our line of sight. 

There are only three known binary pulsar systems in or about our own galaxy 
that will coalesce in less than the age of the universe; of these three, one in 
particular drives the calculation of the rate density. If we attempt to project 
this meager observational data throughout the entire universe, accounting for 
observational biases that cause us to miss some fraction of the actual number of 
binary systems, we find that the rate density of coalescing binaries is 



n 



10~^^^Mpc^^yr~\ (132) 



The great uncertainty in how to project the observations throughout the universe 
is reflected in the factor of 10^ uncertainty in this rate. (For a sense of the 
uncertainties and corresponding controversy surrounding the estimates of the rate 
density of inspiraling binary coalescence, see0B and references therein.) 

The quantity tq, defined in equation |131| , was constructed in such a way that. 



assuming that sources are distributed homogeneously and isotropically throughout 
space, the rate of sources observed above a signal-to-noise p of 8 (p^ of 64) is equal 
to 

N = —rlh. (133) 

For the initial LIGO detectors, tq is only about 13 Mpc; so, the anticipated rate 
of binary inspiral is, even at its most optimistic, low for the first generation of de- 
tectors. Things get better for the more advanced instruments now on the drawing 
board: for these, tq climbs to over 100 Mpc and, by correlating the signal from 
the two LIGO 4 Km detectors and the LIGO 2 Km detector, the effective tq can 
be increased by another factor of 3/2. ^ 

In addition to the SNR, we also observe the scaling of the radiation frequency 



with time. From equation 12S this gives us Ai, the so-called chirp mass, which de- 



pends on the component masses and the source's cosmological redshift. Knowing 



both the SNR and the chirp mass raises the interesting possibihty of measuring 
the Hubble constant: the rate of cosmological expansion. Measuring the SNR 
tells us, up to the complications of the noise and the unknown orientation angles, 
something about the luminosity distance to the source. Similarly, measuring the 
chirp mass tells us, up to the unknown component masses, something about the 
source redshift. There is thus a hidden redshift/luminosity-distance relationship 
in observations of binary inspiral. By statistical analysis of a large number of 
binary inspiral observations, that relationship can be extracted and with it the 
Hubble const ant .ii 



5.1.2 Compact binary coalescence 

Eventually the inspiraling orbit of a binary system with compact components 
must end: the neutron stars collide, the black holes merge, or the one black hole 
tidally disrupts its neutron star companion. The radiation arising from the last 
few inspiral orbits through the coalescence of the compact objects that compose 
the binary may also be a significant source of radiation. 

Unfortunately, very little is really known about the gravitational waves that 
result from the late stages of inspiral or the coalescence of either neutron star 
or black hole binary systems. In both cases the gravitational fields are quite 
strong and dynamical, which would appear to rule-out a perturbative approach 
and require a (numerical) solution to the fully non-linear Einstein field equations.il 
(In the neutron star case, the problem is further complicated by the need to model 
the dynamics of the fluid, which cannot be ignored in a coalescence.) 

Numerical simulations of the coalescence of two black holes in a head-on colli- 
sion have been calculated and provide some guidance: these tend to show that the 
total energy energy is disappointingly small: on order 10~^ of the system's total 
mass-energy.&ii There are also some recent calculations of off-axis collisions, 
which suggest strongly that the maximum energy radiated in such a collision will 
be less than 1% of the total mass energy .H In particular, it is difficult to justify 
the additional, ad hoc factor of 10%Mc^ assumed by some authors in estimating 
the detectability of these sources.. 

Surprisingly, perturbation calculations of the radiation arising from colliding 
black hole spacetimes give results that are in close accord with the limited number 
of fully relativistic numerical simulations that have been performed. 



This accord is difficult to explain and may signal that we have something new to 
learn about the nature of the solutions to the full field equations. 

5.1.3 Black Hole Formation 

Black holes form in the collision of neutron stars at the end-point of neutron 
star binary inspiral; they also form in the core collapse of sufficiently massive 
stars. Unless the formation mechanism is especially symmetric, the new black 
holes that form will be initially quite deformed and will need to radiate away 
their deformations before they can settle down into a quiescent state, which is 
axisymmetric. 

Quiescent black holes are characterized only by their mass M and angular 
momentum J. (And electric charge, too; however, astrophysical black holes are 
unlikely to carry any significant electric charge.) Correspondingly, while the initial 
radiation from the formation of a black hole depends on the details of the forma- 
tion, the final radiation depends principally on M and J. In fact, the late-time 
waveform from a perturbed black hole is a superposition of exponentially damped 
sinusoids, whose frequencies and damping times depend only on M and J, the 
overtone number n and the harmonic order ^ and m of the perturbation. 

Most all of the modes of a black hole are very strongly damped. The most 
weakly damped modes are associated with the fundamental quadrupole-order ex- 
citations. Even these are strongly damped unless the black hole is very rapidly 
rotating. For this reason, we focus attention on the fundamental quadrupole 
modes, which are the most likely to be detectable. Setting aside the start-up 
transient associated with the details of the initial excitation, a good model for the 
"ring-down" of a newly-formed or perturbed black hole is thu^lil 



where the amplitude is averaged (in a root-mean-square sense) over all orientation 
angles. 
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r is the distance from the black hole to the detector, e is the fraction of the total 
mass of the black hole carried away in radiation, and we have assumed that all five 
of the fundamental tone quadrupole modes are excited equally. Corresponding to 
this radiation is an estimated signal-to-noise ratio of 



• an efficiency e for fraction of the rest mass of the system radiated gravita- 
tionally that is equivalent to what is found in black hole collisions, and 

• the effective noise power spectral density is approximately constant over the 
signal bandwidth (which is broad for strongly damped oscillations). 

The rate of black hole formation is entirely uncertain; however, most astro- 
physicists see no reason why the same mechanisms that make neutron stars cannot 
also make black holes at approximately the same rate.0 By our present under- 
standing of formation mechanisms, this rate is not high even at the distance of 
the Virgo cluster (~ 20 Mpc): perhaps as many as a few per year, but likely much 
less. Consequently, should be at least on order 30-35 for a confident detec- 
tion in ideal circumstances.!! The caveat of "ideal circumstances" is an important 
one, however: the character of the waveform for this source — an exponentially 
damped sinusoid — is exactly the kind of technical noise one might expect in a 
real interferometer owing to transient disturbances that affect, for example, the 
suspension of the interferometer mirrors. Thus, without strong assurance that 
what is observed is not a weak disturbance intrinsic to the detector, prospects are 
not good for observing radiation from this source. 

5.1.4 Stellar Core Collapse 

Theoretical models of stellar core collapse, and the corresponding gravitational 
wave luminosity, have a long and checkered history: estimates of the gravita- 
tional wave luminosity have, over the last 30 years, ranged over more than four 
orders of magnitude .1-31-3 It is not simply the luminosity that is unknown: the 
waveforms themselves are also entirely uncertain, leading to a further difficulty 




(139) 



where we have assumed 



in estimating the detectability of this source. (Examples in the hterature can be 
found in the citations.&ll^liiSSiBlli) Nevertheless, it is still possible to evaluate 
what is required of stellar core collapse in order that it be observable in a given 
detector.0 

Suppose that the waveform from supernovae is given by 

h+ = ^«/+m(t) (140) 
r 

= '^f^f^mit) (141) 
r 

a and (3 are constants, /+ and fx are functions of the relative orientation of the 
source with respect to the detector, and m{t) is some function of time which we 
leave undetermined for now. The power radiated into each polarization mode is 
given by 

E+ = «2(/2)M||m|2 (142) 
E. = (3'{fl)Ml\m\\ (143) 

where <> signifies an average over a sphere centered on the source. 

Now assume that equal power is radiated into the two polarization modes. 
Then we can write a and /3 in terms of a single parameter e as 
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In terms of e the power radiated into the + and x polarization states is thus 

Finally, return to consider the time dependence of the waveform m{t). Note 
that, by the Parseval Theorem, 

dt\m{t)\'= fdf{2nff\m{f)\\ (147) 



where m is the Fourier transform of m. Assume, as suggested by numerical simula- 
tions, that the (real) power radiated per unit bandwidth is approximately constant 
for frequencies in the interval (/min, /max) and falls off rapidly outside that band, 
with /min on order 100 Hz and /max approximately 1 KHz. In this approximation, 

(2vr/|m(/)|)^^ . (148) 

/max /mill 



We can now evaluate the SNR we expect from core collapse supernova gravita- 
tional waves incident on the detector. Current calculations suggest e in the range 
1O"^-1O~*M0, with peak power in the 200-300 Hz band.S'SEi For convenience 
here, suppose that the detector noise power spectral density Sh is approximately 
constant over the bandwidth of the signal {i.e., from /min to /max) and (optimisti- 
cally) equal to its approximate value at 100 Hz, and that fmm is very much less 
than /max- The mean-square signal-to-noise ratio is thus 

- 1 , 1 1 
~ 1 H 
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The distance of 15 Mpc is the range to the center of the Virgo Cluster of galaxies: 
if supernovae can be reliably observed to distance, we can expect a rate of approxi- 
mately three per year. Reliable observation of millisecond bursts at a rate of three 
per year, however, requires an SNR of somewhat more than 30. Thus, without 
a very optimistic efficiency e, we can't expect to be able to observe supernovae 
much beyond our own galaxy, and certainly not out to the Virgo cluster. 

5.2 Periodic Sources of Gravitational Radiation 

All the anticipated sources of periodic gravitational waves for the ground-based 
detectors now under construction involve tapping the stored rotational energy of 
rapidly rotating neutron stars. This means that the neutron star must be in some 
way non-asymmetric. The strength of the radiation depends on the degree of 
asymmetry. All the uncertainty associated with periodic sources of gravitational 
radiation arises either with the mechanism for producing the asymmetry or the 
degree of asymmetry. 

At this writing four different kinds of asymmetries, or mechanisms for pro- 
ducing asymmetries, are discussed as possibly leading to detectable gravitational 
waves. We discuss these in subsections ^.2.1| , |5.2.4| , |5.2.5| , |5.2.6| below. Addition- 



ally, we discuss observational constraints on the radiation from isolated pulsars in 



subsection ^.2.2| , and some issues related to the detection of periodic radiation in 



subsection 5.2.3 



5.2.1 Non-axisymmetric rotators. 

All but the youngest neutron stars have a solid, crystalline crust, covered by 
a fluid surface and with a fluid interior. The fluid surface and interior adjust 
themselves to the star's rotation, remaining always axisymmetric and, therefore, 
not contributing to any gravitational radiation. As the star cools or spins-down 
(owing to, e.g., magnetic multipolar radiation if it is a pulsar) the shape of its 
crust cannot adjust continuously to its new conditions. The stresses in the crust 
build until the crust fractures, relieving the stress. The final crust shape is likely 
to be non-axisymmetric and responsible for gravitational radiation as the star 
rotates. 

Suppose that the star is rotating about a principal axis of its moment of inertia 
tensor with rotational rate is /. Let be the moment of inertia along the axis or 
rotation, Ji and I2 be the other two principal moments of inertia, and define e to 
be the difference between Ji and I2 relative to J3: 

e = {h-h)/h. (150) 

Setting aside the very slow spin-down of the system as its angular velocity changes 
and averaging over the angles that describe the relative orientation of the pulsar 
with respect to the detector, the characteristic radiation from this system is given 

by 

hit) ^ /lo cos(47r/t + 0o) (151) 
h„ = -i^'JlIl (152) 

= 4.8 X lO-^^^l^ . (153) 

10^5 g cm-3 10-6 r V300Hzy ^ ' 

The power radiated gravitationally through this mechanism depends, through 
e, on the degree of asymmetry that can be supported by the neutron star crust. 
Alpar and Pines^l have looked at the structure of the crust and the likely strain 
that it can support. For our purposes it is instructive to look at the most ex- 
treme possibility they considered: that the crust is well approximated as a pure 
Coulomb-lattice crust. Such a lattice could sustain a strain some 10^ to 10^ times 
as much as is typical of terrestrial material. When the maximum allowable strain 
is supposed to be supported by the solid part of the neutron star (which is only a 
small fraction of the entire star), one arises with a maximum e of approximately 
10-6. 



This is an extreme value: it depends on the crust being a pure Coloumb lattice, 
assumes that some mechanism has led it to be stressed to its fracture point, and 
that the corresponding strain is principally quadrupolar. For young neutron stars 
one might imagine this conspiracy of circumstance possible; however, for older 
neutron stars plastic flow of the crust would lead to relaxation over the age of 
the most rapidly rotating neutron stars — the so-called millisecond pulsars — , 
reducing the maximum e for these systems to no greater than 4 x 10~^°. 

5.2.2 Observational constraints 

There is observational evidence that, at least for the older, millisecond pulsars, 
e cannot be much greater than this limit. The power L radiated gravitationally 
by a spinning neutron star comes directly from the star's rotation; consequently, 
radiation back-reaction must slow the star in such a way that energy is conserved. 
This leads to a slow spin-down of the star: if P is the spin period, then the rate 
of spin P, assuming that gravitational radiation reaction is the only source of 
spin-down, is 

P = ^. 154 

Since the radiated power L is proportional to jP^ (cf. 0) the measured period 
and period derivative place a strict upper limit on e for isolated pulsars. 

The spin-down rate (P) of most pulsars has been measured. If we take the 
most extreme view and ascribe all of the spin-down to angular momentum carried 
off by gravitational waves, the oblacity of millisecond pulsars still can not exceed 
10~^ for most millisecond pulsars.0 

For a few young, isolated pulsars, timing observations are so good that we can 
place still stronger limits on e: limits that exclude the possibility that a significant 
part of the spin-down is owing to radiation reaction. For these pulsars, not only 
the rate of the spin P but also its second derivative P has been measured. Using 
only that the rotational energy of the star is proportional to J/P^ and that the 
radiated power is proportional to P~" one can quickly show that 

PP 

-^ = \-n. (155) 

If the spin-down is due to quadrupole gravitational radiation reaction, n is equal to 
5 and higher-order radiative moments would lead to larger n. On the other hand. 



if the spin-down is due to, say, magnetic dipole radiation (from the rotation of the 
pulsar's magnetic dipole moment), n is equal to 3. There are no isolated pulsars for 
which the measured PP/P"^ yields an n approaching 5. In particular, for the Crab 
pulsar, which is the prototypical young pulsar, the measured n is approximately 
2.5. This is strongly inconsistent with the suggestion that gravitational radiation 
damping plays an important — let alone dominant — role in the spin-down of the 
Crab, or any other known, pulsar. Thus, we must conclude that effective oblacity 
e of pulsars is, in fact, quite small and the radiation quite weak compared to the 
detector noise. 



5.2.3 Detecting periodic sources 



That the signal amplitude is small is, by itself, not of overwhelming concern. 
Consider for a moment the signal-to-noise associated with an observation of a 
periodic signal, as is given in equation |151| , over a period T long compared to the 
periodicity of the radiation l//gw The SNR is given by 
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Since the signal is monochromatic over the observation period, h is appreciable 
only at frequencies near /gw; consequently, we can set Sh{f) equal to Sh{fgw) 
in equation p.56| . Then, invoking The Parseval Theorem, we can rewrite equation 
1561 in the time domain: 
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i.e., for periodic sources the SNR grows with the observation period. As a result, 
a signal with a small peak amplitude can be detected much more readily if it is 
periodic than if it is a burst signal of finite duration. 

In reality, there are practical limits to how large the SNR can be made. An 
observation can only be so long: for LIGO an observation lasting a full year would 



be quite long. Additionally, the signal is not quite as simple as we have assumed. 
The detector's motion with respect to the source is non-uniform, owing both to 
Earth's rotation about its axis and its orbit about the sun. The time-dependent 
Doppler shift of the signal leads to a frequency modulation, which depends on the 
position of the source in the detector's sky. In order to obtain the growth of 
with time that we found above that frequency modulation must be removed. 

Neither of these issues is significant if we know the position of the source on 
the sky and its frequency. If, on the other hand, we are contemplating a "blind" 
search, across the sky or over a wide bandwidth in frequency, the story changes. 
If we are only interested in one frequency, say /o, then by folding the data we 
can find h{fo). If we are interested in a wide range of frequencies, however, the 
most effective way to find h{f) in that band involves a Fourier transform over 
the interesting bandwidth. For a significant bandwidth (say, several hundred Hz) 
and an observation lasting several months this is a multi-gigapoint FFT. If, in 
addition, we are interested in searching over the sky, then we must demodulate 
differently for different points in the sky. For a several month observation, the 
number of independent patches in the sky is quite large and the computational 
resources required exceed any that might be available now or in the forseeable 
future.Ei Thus, different ways of searching for periodic sources must be developed, 
which will not have the same growth rate of SNR with time. 

Finally — and this is true even if we are interested only in one point on the sky 
and one frequency — as the observation gets longer, we must become concerned 
about the stability of our instrument and the characterization of its noise. One 
cannot expect the noise level to remain stationary over indefinitely long periods; 
additionally, as we observe for longer periods, we are, in an important sense, 
exploring in greater and greater detail the character of the noise in very narrow 
bandwidths. As the observation period increases we must ask, with an increasing 
degree of concern, how confident we are that there are no technical noise sources, 
such as weak, drifting oscillators, that may be masquerading as signal over the 
period of our observation. 

5.2.4 Precession. 

An axisymmetric neutron star, rotating about it symmetry axis, does not radiate 
gravitationally. On the other hand, if the angular momentum is not coincident 



with the symmetry axis — e.g., if the star is precessing — then it will radiate 
gravitationally. Misalignment of an axisymmetric neutron star's angular momen- 
tum and body axes could arise as the result of crustal fractures associated with a 
neutron star quake. 

The same observational constraints that apply to gravitational radiation aris- 
ing from the rotation of a non-axisymmetric neutron star about a principal axis 
also apply to radiation arising from precession of an axisymmetric neutron star 

(cf. §|1D 

While neutron star precession is — in principle — possible, if the neutron star 
is also a pulsar, the precession should also manifest itself as periodic variations 
in the electromagnetic pulse shape. At present there is no observational evidence 
for pulse shape variations induced by free precession. This may be because the 
misalignment is too small to be observed in the pulse shape or because the stresses 
associated with misalignment quickly bring the star back in to alignment. 

Gravitational radiation associated with precession of an axisymmetric star 
occurs at both the rotational frequency and twice the rotational frequencylll; con- 
sequently, it can be distinguished from the radiation associated with a fully non- 
axisymmetric star rotating about a principle axis. Interestingly, observing the 
amplitude of the radiation at both the rotation frequency and twice the rotation 
frequency allows one to determine the all the angles that characterize the orienta- 
tion of the star relative to the line-of-sight (LOS): the angle between the angular 
momentum and the LOS as well as the angle between the body axis and the an- 
gular momentum. If the star is also observable as a pulsar, then one can test 
models of pulsar beaming, since, together with the observed pulse shape, these 
make predictions about the angle between the magnetic axis, the LOS and the 
magnetic axis. 

5.2.5 Thermally driven non-axisymmetry. 

Timing of the x-ray emission from several accreting neutron stars has revealed 
quasi-periodic variability that can be explained as arising from the rapid rotation 
of the underlying neutron star. An intriguing coincidence in these observations 
is that the rotation rate of all these systems appears to be close to equal. This 
suggests that there is some underlying mechanism that insures that accretion 
spins these stars up to — but not beyond — this limiting angular velocity. One 



possibility is that the rotation rate is hmited by gravitational radiation reaction. 

How might gravitational radiation limit the rotation rate of an accreting sys- 
tem? If the accretion leads to a non-axisymmetry in the neutron star then, as the 
star spins-up, the angular momentum radiated by this rotating non-axisymmetry 
increases until it balances the angular momentum accreted, limiting the stars ro- 
tation rate. The angular momentum radiated is, like the radiated power, a strong 
function of angular velocity (J is proportional to Q^, where fl is the angular ro- 
tation rate); so, it is not surprising that the limiting angular velocity should be 
similar for these systems. 

Proposals like this are characteristically made for systems where there appears 
to be some upper limit to the rotation rate. To be plausible, there must be 
some universal mechanism whereby the same process that spins the star up also 
leads to a non-axisymmetry that can cause a radiative loss of angular momentum. 
Recently Bildsten0 offered some promising ideas for a mechanism like this that 
would operate in rapidly accreting, low magnetic field neutron stars like Sco X-1. 

At the core of Bildsten's proposal is the observation that localized heating 
of the neutron star owing to non-isotropic accretion leads to differential electron 
capture rates in the neutron star fluid. These lead, in turn, to density gradients as 
nuclear reactions in the neutron star adjust its composition. Bildsten suggested 
that, if the rotation axis is not aligned with the accretion axis and if some other 
mechanism (in Bildsten's original suggestion, a magnetic field) can break the 
symmetry still further, these density gradients may form in a non-axisymmetric 
fashion. 

There are two big "ifs" in this proposal. Even if the accretion disk is misaligned 
with the star's rotation axis the density perturbations will be distributed symmet- 
rically about the star's rotation axis unless some other mechanism can be shown 
to break the symmetry further. Additionally, though not recognized in the origi- 
nal proposal, buoyancy forces will lead the density perturbations to re-distribute 
themselves in the star symmetrically about its rotation axis, significantly sup- 
pressing the gravitational radiation. For these two reasons the initial excitement 
over the Bildsten proposal has dampened. It should not be extinguished, however: 
the recognition that accretion can lead, through pyro-nuclear reactions, to density 
perturbations that may radiate gravitational is certainly sound and new. With 
time will come greater understanding of where and how this effect may arise in 
nature, and that greater understanding may yet include a robust mechanism for 



producing significant gravitational radiation from accreting neutron stars. 



5.2.6 R-mode instability. 

Also in the last year, Andersson0'Ei discovered a new class of unstable perturbative 
modes of rotating relativistic stars. In the absence of gravitational radiation these 
modes are all stable; however, gravitational radiation back-reaction on the modes 
causes them to undergo exponential growth, feeding off of the rotational energy 
of the star. 

These particular modes are unusual in two different respects: 

1. In the absence of viscosity they are unstable for all angular velocities; 

2. At leading order the radiation is entirely "magnetic" in character: i.e., the 
radiation couples to the momentum density distribution and its time evolu- 
tion, not the mass density distribution and its time evolution. 

Figure ^ gives a schematic view of the character of the fluid velocity perturba- 
tion, relative to the star's rotation, for the lowest order radiating mode (magnetic 
quadrupole) and for low rotational velocity. The lines indicate fluid flow-lines; 
the arrow indicate the relative direction of fluid flow on adjacent flow-lines. At 
low rotation rates there is very little radial component to the fluid motion; so, we 
show only the angular components. 

In the absence of rotation, the fluid in the star can circulate in the pattern 
shown with constant velocity: i. e., there are no restoring forces acting on the fluids 
inertia and the eigenfrequency of the mode is zero. Correspondingly, the magnetic 
quadrupole moment of this fluid mode has constant amplitude. Rotation intro- 
duces Coriolis forces, which act back on the fluid's inertia and cause the circulation 
in these cells to be periodic. The quadrupole moment of the fluid momentum now 
has a second time derivative; correspondingly, it radiates gravitationally. 

The power radiated gravitationally comes, ultimately, from the star's rota- 
tional energy. The radiation thus carries off positive angular momentum from the 
star. The greatest angular momentum can be carried by the modes with azimuthal 
quantum number m equal to ±£; so, focus attention on these modes. Relative to 
the stars rotation, one of these modes is co-rotating and carries positive angular 
momentum and one is counter-rotating and carries negative angular momentum. 
It turns out, however, that both of these modes are dragged, by the stars overall 




Fig. 9. Schematic illustration of the fluid flow-lines for the lowest order r-mode 
perturbation relative to the star's rotation. The arrow indicate the relative di- 
rection of fluid flow on adjacent flow-lines. In the upper panel the stars surface 
(co-latitude 9, longitude 0) is mapped to the plane; in the lower panel, the view 
is looking down along the stars rotation axis. In both cases, the flow pattern is 
dragged along by the star's rotation. At low rotation there is essentially no radial 
component to the fluids motion; so, we show only the angular components. 



rotation, so that viewed from an inertial frame both carry positive angular mo- 
mentum. The back-reaction of the radiation adds negative angular momentum 
to each of these modes; correspondingly, the co-rotating mode is damped and 
the counter- rotating mode is anti-damped: i.e., it grows. As it grows, of course, 
it radiates more strongly, leading to greater anti-damping: the mode undergoes 
exponential growth. 

This general mechanism by which gravitational radiation reaction can lead 
to amplification of a mode that is counter-rotating in the body frame but co- 
rotating in the inertial frame was first discussed by Zel'dovich in the context of 
rotating black holes. His concern was that, through this mechanism, a rotating 
black hole should radiate — a result that anticipated the more general result of 
Hawking. The first application to stars came from Chandrasekhar@ and Friedman 
and SchutJi whose focus, however, was on a different set of modes. 

Viscosity damps the R-modes. If the viscosity is large enough, then the viscous 
damping exceeds the anti-damping caused by radiation reaction and the mode is 
stabilized. Present understanding of the viscosity of neutron star fluid suggests 
that there is a short period in the life of a new-born neutron star, lasting perhaps 
one year, when the mode may be unstable and a rapidly rotating star may be a 
strong, nearly periodic source of gravitational radiation. The radiation is "nearly" 
periodic because its amplitude is so great that, over the course of the year, the 
stars angular rotation rate may evolve from 10^ to 10^ Hz, simply due to the 
angular momentum carried away in the radiation. 



5.3 Stochastic Gravitational Radiation 

In the previous subsections we discussed burst and periodic sources of gravitational 
radiation. In both cases the discussion focused on the source and the radiation 
was characterized typically in terms of a waveform h{t)., which depends on the 
details of a source and its orientation with respect to the detector. 

The situation for a stochastic gravitational wave signal is different. A stochas- 
tic gravitational wave signal is intrinsically random in character. In particular, it 
is not generated by an isolated source, it is not incident on the detector from a 
single direction, and it has no characteristic waveform. 

In fact, a stochastic signal can be treated as just another source of detector 
noise. The stochastic radiation has an h{t) that is characterized solely by its 



correlation functions or associated spectra; correspondingly, the detector output 
owing to the stochastic signal is characterized in terms of correlation functions or 
associated spectra. 

Detection of any signal hinges on observing some characteristic that distin- 
guishes signal from noise. If a stochastic signal appears in a detector no different 
than intrinsic detector noise, how do we make the critical distinction that allows 
us to say we have detected a signal? 

The essential difference between the action of a stochastic signal and intrinsic 
detector noise is that stochastic radiation incident on two detectors is correlated, 
and the correlation depends in a completely predictable way on the relative orien- 
tation and separation of the detectors. Any gravitational wave signal — stochastic 
or otherwise — can be resolved into a superposition of plane waves of different 
frequencies and propagation directions bathing the detectors. A component of the 
radiation of given frequency and incident direction drives two or more detectors 
coherently, with a phase delay that depends on the incidence direction, detec- 
tor separation and radiation wavelength. For components whose wavelength is 
much larger than the separation between the detector, the phase difference is only 
weakly dependent on the radiation wavelength or incident direction; so, summed 
over incident direction there is a strong correlation between the output of the 
detectors. On the other hand, for radiation components with wavelengths much 
smaller than the separation between the detectors the phase difference depends 
strongly on the incident direction and the wavelength; so, summed over incident 
direction, the correlation between the output of the detectors is weak. Thus, in 
the presence of a stochastic gravitational wave background the output of two or 
more detectors should show predictable correlations that depend on their relative 
separations, relative orientations, and the stochastic signal's spectrum.il 

When considering sources of detectable stochastic gravitational radiation for 
ground-based detectors, it is conventional to enumerate the contributions of pri- 
mordial origin: e.g., radiation arising during an inflationary epoch in the early 
universe,ii from the decay of a cosmic string network]^ or from a phase transition 
in the early universe.&il Less frequently discussed — perhaps because it is so 
mundane — are the contributions arising from the confusion limit of discrete but 
unresolved sources: e.g., core-collapse supernovaeS or binary inspiral. In fact, 
the contribution at low-frequencies (periods of hours to minutes) to the stochas- 
tic signal from unresolved galactic binary systems is expected to be many times 
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# events 


"rate" 




0-5 s 


3 


0.6/s 


number m 


5-10 s 


4 


0.6/s 


any interval 


10-15 s 


11 


0.6/s 


varies from 


15-20 s 


2 


0.6/s 


long-run rate 


20-25 s 


5 


0.6/s 




0-25 s 


25 


1.0/s 


long-run rate 



Table 1. Number of Poisson distributed events, with mean rate 1/s, occurring in 
consecutive five second intervals. Note that the number of events in any given 
interval varies: the rate is only the mean of the distribution of "rates" calculated 
over many intervals. 

greater than the intrinsic detector noise of the proposed space-based interferomet- 
ric detector LISA.iiil 

Let's talk a bit more about how discrete, well-defined sources — whose wave- 
forms may be well known — can superpose to form a stochastic signal of pre- 
dictable spectrum. (We focus on burst sources here; however, a similar story can 
be told for periodic ones.il) 

Astrophysical sources of gravitational-wave bursts strong enough to be ob- 
served by conceivable detectors are certainly independent events. Setting aside 
the slow evolution of the source population on cosmological timescales, the num- 
ber of events whose radiation impinges upon the receiver during any finite length 
observation is Poisson distributed: if the event rate is and the observation 
period is r, then 



P{n\TN) = 




(159) 



It is critical to recognize that the number of events exciting the detector at 
any given moment is not the constant equal to the product of the event rate and 
the signal duration. The actual number of events in the observation period varies. 
Averaged over many intervals it has a mean value, which is the product of the 
event rate and the signal duration. For example, table |l| shows the number of Pois- 
son distributed events over consecutive five second intervals, the rate calculated 
separately in each interval, and the overall average rate. 

Now let's suppose that the "events" we are discussing are the arrival at our 
detector of the initial wavefront associated with an impulsive perturbation of a 
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Fig. 10. The prototypical waveform associated with an impulsive perturbation of 
a black hole: an exponentially damped sinusoid. 



black hole. That waveform is a damped sinusoid: 

h{t) oc e-2-/(*-*o)/Q sin 27r/(t - to) for t > to. 



(160) 



The signal from a single, prototypical event is shown in figure |T0|. The left-hand 



panels of figure |Tl], on the other hand, shows h(t) at the detector when signals 
exactly like these arrive at the detector with different rates. In the top panel the 
signals are, for the most part, clearly distinguishable. In the middle panel the rate 
is higher and separate signals can be distinguished only in exceptional cases. In 
the bottom panel the rate is higher still and the identity of the individual signals 
is completely obscured. 

Especially in the bottom panel of figure 0, the signal h{t) at any given t is the 
superposition of a large, but random, number of signals. Similarly, the amplitude 
of each contributing signal is itself random (corresponding to the signal amplitude 
at a random moment relative to the signal start time). In circumstances like these 
we expect The Central Limit Theorem^ to apply, leading to a normal distribution 



for h(t). On the other hand, at the lower rate pictured in the top panel, there 
are not enough events superposed at any given moment for us to expect h{t) to 



exhibit a normal distribution. In the right-hand panels of figure |TT] we show the 
distribution of h{t) taken from the left-hand panel and find that these intuitions 
are borne out. 

Note that none of these conclusions depend in any way on the details of the 
signal from an individual source: instead of the superposition of damped sinu- 
soids we could just as well have constructed h{t) from the superposition of binary 



inspiral signals in a fixed bandwidth. We do exactly that in figure 0. Note how 
the distributions of h{t) at high rate are, for superpositions of large numbers of 
signals, identical {i.e., they are normal distributions). 

Despite the fact that the distribution of h{t) arising from the superposition of a 
large number of sources are identical, visual inspection suggests that there are still 
differences. These differences are associated with the correlations: the distribution 
of the products h{t)h{t + r) as a function of r. These are very different for the 
damped sinusoids, which are characteristic of black hole perturbations, and the 
"chirps" , which are characteristic of binary inspiral. In figure O we show the 



power spectral densities of the superpositions in figures |TI| and |T^. Note how 
the power spectral densities of the stochastic signals formed from the random 
superposition of events of a given character preserve the spectral shape of the 
underlying signal, with the overall amplitude proportional to the event rate. It 
is this property of the superposition that permits us to predict the character 
of the stochastic signal arising from the confusion limit of a large number of 
sources; conversely, observation of a stochastic signal provides us, through its 
amplitude and spectrum, information about the underlying source, its rate and 
spatial distribution. 

This last point gives us the prospect of using a detected stochastic background, 
arising from the superposition of unresolved sources, to perform a source "cen- 
sus": a determination of the density of sources in space. This information, in 
turn, can shed light on astrophysical process that are otherwise unobservable. An 
excellent example of this comes from LISA observations of the background from 
close white dwarf binary systems (CWDBs).iil These binary systems, which have 
orbital periods ranging from days to hours, are so close that they are optically 
unobservable as binaries. We know they exist, nevertheless, because we can see 
their progenitor systems. CWDBs arise as one of the end-points of binary star 




Fig. 11. The left-hand panes show waveforms corresponding to the superposition 
of many discrete and ideahzed black hole formation events, with the number of 
events in a fixed interval Poisson distributed. The rate increases by an order of 
magnitde from the top to the middle panel, and again from the middle to the 
bottom panel. The right-hand panes show the distribution of wave amplitude 
derived from the corresponding left-hand panel. 
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Fig. 12. As in figure 
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, but with band-limited binary inspiral waveforms. 



evolution, following a particularly difficult to understand and model evolutionary 
stage where the two stars are orbiting each other within a single envelope of gas. 
Once they emerge from this final stage of common envelope evolution they orbit 
each other as if point masses. Gravitational radiation reaction leads to a slow, 
secular decay in the orbit of these systems until they become so close that mutual 
tidal interactions lead to rapid orbital decay or stellar disruption. 

Let dn/df be the number density of CWDBs with orbital frequency /. The 
amplitude of the stochastic gravitational wave signal from this population of ob- 
jects is, at frequency 2/, proportional to dn/df. CWDBs are introduced into 
the population at a rate dn+/df as they emerge from the final stage of common 
envelope evolution at orbital frequency /; similarly, they leave the population at 
a rate dn_/df as they disrupt at orbital frequency /. Knowing the rate dn^/df 
would provide us valuable information about common envelope evolution, which 
we cannot obtain through optical observations; similarly, the rate dn^/df depends 
on the white dwarf mass spectrum (we know the white dwarf equation of state 
quite well) and knowing it would reveal the mass spectrum of binaries emerging 
from the common envelope evolution stage. 

We can determine both dn^/df and dn^/df from LISA observation of the 
CWDB stochastic signal spectrum. Observation of this spectrum determines, 
as we have seen, dn/df. Once injected into the population, a CWDB's orbit 
evolves owing exclusively to gravitational radiation reaction, which proceeds at 
the rate df/dt, until it is removed from the population through coalescence or 
disruption. In steady state we thus have a continuity equation governing the 
CWDB population: 



The orbital frequency evolution rate df/dt is known (it is proportional to /^^^'^) 
and LISA observations will determine the number density dn/df; consequently, 
from LISA observations we can determine dfij^/df and dn^/df and learn about 
the end of common envelope evolution and the mass spectrum of white dwarfs in 




(161) 
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Fig. 13. Power spectra corresponding to the superpositions of black hole ringdown 
signals in figure |TT] and binary inspiral signals in figure |T2|. The power spectra 
reflects the underlying source of the stochastic signal (black hole ringdown or bi- 
nary inspiral) while the amplitude is proportional to the rate of individual sources 
contributing to the signal. 



6 Conclusions 



In these lectures I've tried to give a brief overview of how we think about gravi- 
tational waves when we set out to detect them, and provide a snapshot of current 
thinking on the anticipated wave sources. Along the way, I've tried to describe 
some of the science we can hope to do once we can reliably detect gravitational 
wave sources. 

The principal difficulty in discussing the sources that we hope to observe is 
our real lack of knowledge of their character. As is often the case, however, this 
difficulty is really a disguised opportunity: when the detectors come on-line and we 
begin to detect gravitational radiation sources, we will not simply be confirming 
what we already know, but learning things entirely new about the cosmos! 
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